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Abstract 

The  goals  of  this  thesis  were  to  quantify  the  behavior  of  gastropod  larvae  (mud  snails 
7 Jyanassa  obsoleta)  in  turbulence,  and  to  investigate  how  that  behavior  affects  larval  supply 
in  a  turbulent  coastal  inlet.  Gastropod  larvae  retract  their  velums  and  sink  rapidly  in 
strong  turbulence.  Turbulence- induced  sinking  would  be  an  adaptive  behavior  if  it  resulted 
in  increased  larval  supply  and  enhanced  settlement  in  suitable  coastal  habitats. 

In  laboratory  experiments,  mud  snail  larvae  were  found  to  have  three  behavioral  modes: 
swimming,  hovering,  and  sinking.  The  proportion  of  sinking  larvae  increased  exponentially 
with  the  turbulence  dissipation  rate  over  a  range  comparable  to  turbulence  in  a  tidal  in¬ 
let,  and  the  mean  larval  vertical  velocity  shifted  from  upward  to  downward  in  turbulence 
resembling  energetic  nearshore  areas. 

The  larval  response  to  turbulence  was  incorporated  in  a  vertical  advection-diffusion 
model  to  characterize  the  effects  of  this  behavior  on  larval  supply  and  settlement  in  a  tidal 
channel.  Compared  to  passive  larvae,  larvae  that  sink  in  turbulence  have  higher  near-bed 
concentrations  throughout  flood  and  ebb  tides.  This  high  larval  supply  enables  behaving 
larvae  to  settle  more  successfully  than  passive  larvae  in  strong  currents  characteristic  of 
turbulent  tidal  inlets. 

A  study  was  conducted  at  Barnstable  Harbor,  MA  to  estimate  the  responses  of  lar¬ 
vae  to  turbulence  in  the  field.  Gastropod  larvae  from  different  coastal  environments  had 
genus-specific  responses  to  turbulence,  suggesting  that  larvae  use  turbulence  for  large-scale 
habitat  selection.  On  ebb  tides,  mud  snail  larvae  had  a  similar  response  to  turbulence  as 
in  the  laboratory,  with  greater  sinking  velocities  in  strong  turbulence.  Behavior  estimates 
differed  for  flood  and  ebb  tides,  indicating  that  additional  physical  cues  influence  behavior. 
Turbulence-induced  sinking  behavior  would  enhance  retention  and  promote  settlement  of 
mud  snail  larvae  in  habitats  like  Barnstable  Harbor. 

Thesis  Supervisor:  Lauren  S.  Mullineaux 

Title:  Senior  Scientist 

Woods  Hole  Oceanographic  Institution 
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Chapter  1 


Introduction 


Many  benthic  invertebrates  disperse  via  planktonic  larvae  whose  behavior  potentially  in¬ 
fluences  large-scale  settlement  patterns.  Larval  dispersal  and  settlement  are  the  result  of 
coupled  physical  and  biological  processes.  Larval  behavior  (swimming  or  sinking)  allows 
larvae  to  transport  themselves  vertically,  and  vertical  movement  can  affect  the  horizontal 
transport  of  larvae  by  positioning  them  in  particular  currents.  Vertical  swimming  or  sinking 
also  affects  the  supply  of  competent  larvae  to  bottom  substrates  prior  to  settlement.  Larvae 
that  respond  opportunistically  to  environmental  cues  are  more  likely  to  settle  in  suitable 
coastal  habitats. 

Gastropod  larvae  sink  in  turbulence,  and  this  behavior  is  expected  to  affect  the  supply 
of  larvae  to  benthic  habitats.  I  define  larval  supply  as  the  concentration  of  larvae  near  the 
bottom,  because  near-bottom  larvae  are  available  for  settlement.  Larval  supply  is  influenced 
not  only  by  behavioral  responses  to  small-scale  turbulence,  but  also  by  larger  scale  vertical 
mixing.  These  complex  biophysical  interactions  control  the  delivery  of  larvae  to  benthic 
habitats.  My  goal  was  to  investigate  the  interactions  between  turbulence,  larval  behavior, 
and  larval  supply  of  intertidal  gastropods. 

1.0.1  Mud  Snails 

Mud  snails  ( Uyanassa  obsoleta)  were  used  as  a  model  organism  because  they  are  significant 
and  conspicuously  abundant  (Fig.  1-1)  members  of  the  intertidal  community.  Mud  snails 
are  benthic  engineers  -  they  rework  sediments,  ingest  recently-settled  larvae  [4],  and  alter 
the  distributions  of  other  organisms  through  disturbance  [2,  3]  and  competition.  Mud 
snails  are  also  infamous  as  the  intermediate  hosts  of  many  duck  and  fish  parasites  [1],  and 
as  unwanted  invaders  of  west  coast  habitats  [5]. 
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Figure  1-1:  (a)  Mud  snail  aggregations  on  intertidal  flats  of  Barnstable  Harbor,  MA  and 
(b)  mud  snail  larvae  that  are  ready  to  settle. 

Mud  snail  larvae  spend  1.5  to  4  weeks  in  the  plankton  before  settling.  When  mud  snail 
larvae  are  competent  to  metamorphose,  they  can  select  substrates  for  settlement  over  small 
scales  [6],  Behaviors  of  mud  snail  larvae  in  the  water  column  were  unknown.  This  thesis 
describes  the  behavior  of  mud  snail  larvae  in  turbulence,  and  the  potential  effects  of  this 
behavior  on  larval  supply  and  settlement. 

1.0.2  Turbulence 

Two  aspects  of  turbulence  are  important  in  this  thesis.  The  first  is  turbulence  at  the 
larval  scale  (mm7s).  Larvae  interact  with  the  smallest  eddies,  described  by  the  Kolmogorov 
microscales  rj  (length  scale),  r  (time  scale),  and  v  (velocity  scale).  These  scales  are  related 
to  the  kinematic  viscosity  v  and  the  energy  dissipation  rate  e  [7]  by 


/j/3y 

(1.1) 

-  -  & 

(1.2) 

V  =  (i/s)4  . 

(1.3) 

The  dissipation  rate  quantifies  the  transfer  of  turbulent  kinetic  energy  to  heat  energy. 
Although  energy  is  generated  at  large  scales,  it  is  passed  to  smaller  scales  when  large  eddies 
shed  smaller  eddies  and  so  on  (the  energy  cascade  [e.g.,  7]).  The  smallest  eddies  contain  too 
little  energy  to  overcome  viscous  forces,  and  the  remaining  energy  is  dissipated  as  heat.  The 
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dissipation  rate  is  related  to  the  total  energy  produced  as  well  as  the  scales  of  the  smallest 
eddies.  I  use  e  throughout  this  thesis  to  quantify  the  turbulence  experienced  by  larvae. 

The  second  important  aspect  of  turbulence  is  the  vertical  mixing  of  larvae  through  the 
water  column.  Most  larvae  are  able  to  swim  or  sink,  but  all  larvae  are  diffused  by  turbulent 
eddies.  Turbulent  mixing  of  particles  by  eddies  is  described  by  the  eddy  diffusivity  K .  The 
relative  strength  of  diffusivity  vs.  behavioral  advection  (swimming  or  sinking)  determines 
the  extent  that  behavior  influences  larval  distributions. 

1.0.3  Thesis  Overview 

The  goals  of  this  thesis  were  to  quantify  the  sinking  response  of  mud  snail  larvae  as  a  function 
of  turbulence,  and  to  estimate  the  effects  of  this  behavior  on  larval  supply.  In  Chapter  2,  I 
use  a  mixture  model  to  quantify  the  behavior  of  competent  larvae  in  laboratory-generated 
turbulence.  In  Chapter  3,  I  use  an  advection-diffusion  model  to  characterize  the  effects  of 
turbulence-induced  sinking  behavior  on  larval  supply  and  settlement  in  tidal  channels.  In 
Chapter  4,  I  fit  an  advection-diffusion  model  to  larval  distributions  in  a  tidal  channel  to 
estimate  larval  responses  to  turbulence  in  the  field. 
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Sinking  behavior  of  gastropod  larvae  ( Ilyanassa  obsoleta)  in  turbulence 

Heidi  L  Fuchs,1  Lauren  S.  Mullineaux,  and  Andrew  R.  Solow 

Woods  Hole  Oceanographic  Institution,  Woods  Hole,  Massachusetts  02543 

Abstract 

Larvae  of  coastal  gastropods  sink  in  turbulence  and  may  use  nearshore  turbulence  as  an  initial  settlement  cue 
Our  objective  was  to  quantify  the  relationship  between  turbulence  and  the  proportion  of  sinking  larvae  for  competent 
mud  snail  veligers  ( Ilyanassa  obsoleta).  We  exposed  larvae  to  a  range  of  field-relevant  turbulence  conditions  (e  = 
8.1  X  10~3  to  2.7  X  10°  cm2  s-3)  in  a  grid-stirred  tank,  holding  other  factors  constant.  We  used  a  video  plankton 
recorder  to  record  larval  movements  in  still  water  and  in  turbulence.  Larval  trajectories  and  velocity  measurements 
were  extracted  using  video-image  analysis.  We  also  measured  turbulent  flow  velocities  independently,  using  laser 
Doppler  velocimetry.  To  interpret  empirical  measurements  in  terms  of  larval  behavior,  we  developed  a  three- 
component,  normal  mixture  model  for  vertical  velocity  distributions  of  larvae  in  turbulence.  The  model  was  fitted 
to  observed  larval  velocities  by  maximum  likelihood,  to  estimate  the  proportions  of  sinking,  hovering,  and  swim¬ 
ming  larvae.  Over  the  range  of  turbulence  intensities  found  in  typical  coastal  habitats,  the  proportion  of  sinking 
larvae  increased  expooentially  (r2  =  0.89)  with  the  log  of  the  turbulence  dissipation  rate.  The  net  mean  behavioral 
velocity  of  the  larvae  shifted  from  positive  to  negative  when  the  dissipation  rate  reached  — 10“ 1  cm2  s~3.  By  sulking 
when  they  enter  turbulent,  shallow  water,  competent  larvae  could  improve  their  chances  of  settling  in  favorable 
coastal  habitats. 


Very  little  is  known  about  how  larval  behavior  in  the 
plankton  affects  patterns  of  larval  supply  and  settlement  of 
benthic  invertebrates.  Much  work  has  been  done  to  describe 
larval  behavior  during  the  exploration  of  substrates,  when 
larvae  can  sometimes  select  settlement  sites  over  small 
scales  (millimeters  to  centimeters).  Less  progress  has  been 
made  on  understanding  the  behavioral  contribution  while  lar¬ 
vae  are  transported  through  the  water  column  to  benthic  hab¬ 
itats.  Under  some  hydrodynamic  conditions,  larvae  could 
settle  more  successfully  if  they  responded  to  waterborne  cues 
by  sinking  toward  the  benthos.  If  the  swimming  velocity  and 
gravitational  sinking  velocity  differ  by  a  factor  of  two  or 
more,  behavioral  changes  can  significantly  affect  larval  sink- 
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ing  fluxes  (Gross  et  al.  1992;  Eckman  et  al.  1994)  and  the 
time  it  takes  for  larvae  to  reach  the  bottom  (McNair  et  al. 
1997).  Yet,  with  a  few  exceptions  (e  g.,  Pawlik  and  Butman 
1993;  Tamburri  et  al.  1996;  Welch  and  Forward  2001),  it  is 
unknown  whether  larvae  change  their  behavior  in  response 
to  conditions  in  the  water  column. 

Ciliated  larvae  are  generally  assumed  to  reach  the  bottom 
boundary  layer  as  passive  particles  (Butman  1987;  Abclson 
and  Denny  1997),  but  passive  deposition  alone  cannot  ex¬ 
plain  some  population  distributions.  In  Barnstable  Harbor, 
Massachusetts,  mud  snails  ( Ilyanassa  obsoleta )  are  the  most 
conspicuously  abundant  megafauna  on  the  intertidal  mud 
flats,  yet  their  settlement  in  the  harbor  seems  improbable, 
because  swimming  larvae  have  mean  upward  velocities  (H. 
L.  Fuchs  unpubl.  data).  Although  swimming  larvae  are  un¬ 
likely  to  be  deposited  on  the  bottom,  sinking  larvae  have  a 
good  chance  of  sealing  along  with  fine  sand  in  the  harbor, 
and  a  behavioral  switch  may  explain  the  apparently  suc¬ 
cessful  settlement  of  mud  snail  larvae. 

With  heavy  shells  for  ballast,  gastropod  larvae  can  alter 
their  vertical  flux  by  changing  their  mode  of  behavior  from 
swimming  to  sinking.  Veligers  are  weak,  ciliary  swimmers 
with  dense  shells,  and  their  gravitational  sinking  velocities 
are  greater  in  magnitude  than  their  maximum  swimming  ve¬ 
locities  in  any  direction  (e.g.,  Hidu  and  Haskin  1978).  When 
larvae  switch  from  swimming  to  sinking  mode,  the  advective 
component  of  vertical  flux  changes  accordingly.  Most  settle¬ 
ment  models  assume  a  constant  larval  velocity,  ignoring  pos¬ 
sible  effects  of  behavioral  changes  on  larval  fluxes  (but  see 
Eckman  et  al.  1994).  Yet  laboratory  observations  (Crisp 
1984;  Young  1995)  and  field  evidence  (Barile  et  al.  1994) 
suggest  that  gastropod  larvae  pull  in  their  velums  and  sink 
when  disturbed,  and  this  behavioral  change  may  affect  set¬ 
tlement  dynamics. 

Turbulence  could  provide  an  initial  cue  for  larvae  to  sink 
and  explore  for  settlement  sites  (Chia  et  al.  1981).  Many 
larvae  settle  preferentially  on  particular  sediments  (Snel- 
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Table  1.  Re  presen  tali  ve  dissipation  rates  for  ocean  regions. 


Location 

e  (cm-  s  *) 

Source 

Open  ocean  (mixed  layer) 

10  MO*2 

Dillon  and  Caldwell  (1980) 

Continental  shelf  (mixed  layer) 

10-M0~2 

Oakey  and  Elliott  (1982) 

Tidal  channels  and  estuaries 

10  M0° 

Gross  and  Nowell  (1985) 

Surf  z,one 

10-M02 

George  et  al.  (1994) 

grove  et  al  1998)  or  in  the  presence  of  conspecifics  (Schel- 
tema  et  al.  1981)  or  chemical  cues  (Pawlik  1992;  Tamburri 
et  al.  1996).  However,  preferred  settlement  sites  can  only  be 
detected  near  the  bottom  and  over  small  spatial  scales  (tens 
of  centimeters).  In  shallow  nearshore  areas,  sinking  larvae 
would  have  more  contact  with  the  bottom  and  more  settle¬ 
ment  opportunities  than  swimming  larvae.  There  would  be 
no  settlement-related  benefits  for  larvae  that  sink  in  deeper, 
offshore  areas,  away  from  suitable  habitats.  Turbulence  dis¬ 
sipation  rates  generally  increase  from  offshore  to  inshore 
regions  (Table  1)  and  might  indicate  to  larvae  when  they  are 
entering  potential  habitat  areas.  We  hypothesize  that  turbu¬ 
lence  above  some  threshold  level  provides  a  primary  settle¬ 
ment  cue  for  mud  snail  larvae  and  that  larvae  respond  to 
this  cue  by  sinking  more  frequently.  This  is  a  behaviorally 
mediated  deposition  hypothesis:  larvae  are  deposited  in  near¬ 
shore  environments  because  the  behavioral  response  to  a  hy¬ 
drodynamic  cue  increases  their  sinking  fluxes  in  coastal  ar¬ 
eas. 

We  conducted  larval  behavioral  experiments  in  a  grid- 
srinred  turbulence  tank  of  the  type  used  extensively  in  re¬ 
search  on  fluid  turbulence  (e.g.,  Hopfinger  and  Toly  1976; 
De  Silva  and  Fernando  1994)  and  plankton  feeding  rates  in 
turbulence  (e.g.,  Saiz  1994;  MacKenzie  and  Kiprboe  1995). 
Grid  tanks  are  ideal  for  understanding  larval  behavior  in 
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Fig.  1.  Schematic  of  the  rurbulence  tank  as  viewed  from  the  side. 
The  small  rectangle  indicates  the  location  of  video  measurements, 
centered  19.1  cm  from  the  bottom  and  11  cm  from  the  outer  walls. 
Asterisks  indicate  the  locations  of  LDV  measurements. 


water-column  turbulence,  because  they  lack  a  developed 
boundary  layer  and  substrate-related  settlement  cues.  Larvae 
are  not  induced  to  explore  or  attach  to  substrates  in  such 
tanks,  and  their  observed  activities  are  representative  of  be¬ 
havior  in  the  water  column. 

Materials  and  methods 

Larval  cultures — Mud  snail  (/.  obsoleta)  egg  capsules 
were  collected  at  Barnstable  Harbor  on  30  June  2002  and 
divided  into  12-liter  buckets  of  filtered  seawater  in  a  20°C 
culture  room.  Larvae  hatched  out  over  10  d  and  were  sieved 
daily  into  fresh  culture  buckets.  Cultures  were  continually 
aerated,  and  seawater  was  changed  every  other  day.  The  lar¬ 
vae  were  fed  —lO5  cells  ml"1  of  lsochrysis  galbana  and 
Thalassiosira  pseudonana .  Experiments  were  done  in  the 
culture  room,  so  that  larvae  experienced  consistent  environ¬ 
mental  conditions.  Experiments  began  when  larvae  reached 
24  d  of  age  and  were  competent  to  metamorphose. 

Turbulence  tank—' Our  experiments  were  done  in  a  1 03- 
liter,  grid-stirred  turbulence  tank  (Fig.  1).  The  grid  was  cen¬ 
tered  at  40  cm  from  the  bottom  of  the  tank  and  1 0  cm  from 
the  free  surface,  and  was  stirred  from  below  with  an  oscil¬ 
lation  amplitude  of  11.34  cm.  McKenna  (2000)  provides  a 
detailed  description  of  the  tank.  All  measurements  were 
made  at  a  point  far  enough  from  the  grid  to  be  in  an  area 
of  homogeneous,  nearly  isotropic  turbulence  (De  Silva  and 
Fernando  1994),  and  as  far  as  possible  from  any  boundaries. 
Turbulence  intensity  was  proportional  to  grid  oscillation  fre¬ 
quency,  which  was  controlled  by  setting  the  voltage.  The 
turbulent  Reynolds  number  Rem  ( see  Table  2  for  a  descrip¬ 
tion  of  symbols )  was  calculated  as  in  Hopfinger  and  Toly 
(1976),  with  empirical  constants  given  in  McKenna  (2000). 
Treatment  levels  were  Rem  =  (100,  200,  300,  400,  500, 
600),  corresponding  to  the  turbulence  dissipation  range  ex¬ 
pected  in  tidal  inlets.  In  our  measurement  volume,  there  was 
a  slight  (<10  mm  s~l)  downward  mean  flow  due  to  the  weak 
circulation  generated  by  grid  motion  (McDougall  1979;  Mc¬ 
Kenna  2000).  Mean  flow  was  accounted  for  in  our  analysis 
of  larval  behavior  (see  “Analysis”  section). 

Experimental  design — Behavioral  experiments  were  rep¬ 
licated  six  times  with  24-d-old,  competent  larvae,  and  no 
larvae  were  reused.  For  the  last  replicate,  only  27-d-old  lar¬ 
vae  were  available.  In  each  replicate,  the  turbulence  tank  was 
filled  with  0.2  pan  filtered  seawater,  at  a  temperature  within 
±1°C  of  room  temperature.  Several  thousand  larvae  were 
gently  added  to  the  tank,  along  with  —2.5  X  104  cells  ml-1 
of  food.  Because  of  human  error,  no  food  was  added  in  rep¬ 
licate  3. 
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Table  2.  List  of  symbols.  All  velocities  are  vertical,  with  positive 
upwards. 


Symbol  Description 

l  Eddy  length  scale 

p,  Probability  density  of  wL  in  mode  i  in  still  water 

P  Total  probability  density  of  wL  in  still  water 

pr  Probability  density  of  vvL  in  mode  i  in  turbulence 

Pr  Total  probability  density  of  wL  in  turbulence 

Re^jy  Theoretical  turbulent  Reynolds  number 

W9  Flow  velocity  [(~A(iu,0,  ff^)] 

W(  Larval  behavioral  velocity  in  mode  i  of)] 

Wr  Larval  relative  velocity  in  mode  i  [(~/V(/xr,  of,)] 
w  Observed  flow  velocity 

(w)  Mean  flow  velocity  (=/i,) 

wL  Observed  larval  velocity 

(wL)  Net  mean  behavioral  velocity 

Behavioral  velocity  range 
w'  RMS  flow  velocity  from  direct  estimate 

vv'  RMS  flow  velocity  from  spectral  estimate  (*sif<1) 

RMS  flow  velocity  from  empirical  relationships 
RMS  flow  velocity  from  maximum  likelihood  {  =  a0) 
a,  Proportion  of  larvae  in  mode  i  in  still  water 

ctyi  Proportion  of  larvae  in  mode  /  in  turbulence 

e  Turbulence  dissipation  rate 

tj  Kolmogorov  length  scale 

v  Kinematic  viscosity  of  seawater  (=0  01  cm-  s~‘) 

of  Measurement  noise  variance 


In  each  replicate,  the  larvae  were  exposed  to  alternating 
periods  of  calm  water  and  turbulence  (Fig.  2).  After  an  initial 
acclimation  period,  video  measurements  were  collected  with 
no  flow  in  the  tank.  Then  six  turbulence  treatments  were 
administered  in  a  randomized  sequence  determined  by  a  Lat¬ 
in  square.  Rest  periods  between  treatments  allowed  larvae 
to  regain  their  calm-water  behavior.  Video  measurements 
were  collected  using  an  analog  video  plankton  recorder.  The 
video  frame  was  3X4  cm  with  an  estimated  2-cm  depth  of 
field  and  was  illuminated  by  an  infrared  spotlight  (A  =  800 
nm).  Larvae  were  already  acclimated  to  the  ambient  room 
lighting,  so  no  light-related  behavioral  changes  were  ex¬ 
pected  during  the  experiments. 

Before  each  replicate,  a  subsample  of  larvae  was  removed 
from  the  culture.  Larval  competency  was  verified  by  putting 
larvae  in  a  petri  dish  with  Barnstable  Harbor  sediment 
(Scheltema  1961);  most  slopped  swimming  and  burrowed 
into  the  sediment  within  3.5  h  and  metamorphosed  within 
24  h  (n  =  25).  Other  larvae  were  killed  with  a  few  drops  of 
ethanol  for  fall-velocity  measurements.  Dead  larvae  with  re¬ 
tracted  velums  were  pipetted  gently  into  the  top  of  a  2-liter 
glass  cylinder  (45  cm  tall,  7.7  cm  diameter)  filled  with  room- 
temperature  seawater.  Larvae  were  timed  ( n  =  20)  as  they 
fell  through  a  12.5-cm  distance  near  the  bottom  of  the  cyl¬ 
inder.  We  also  anaesthetized  larvae  using  MgCl2  and  at¬ 
tempted  to  measure  their  fall  velocities  with  velums  extend¬ 
ed,  but  these  larvae  often  recovered  their  mobility  in  the 
seawater  column.  Larval  shell  lengths  were  measured  {n  = 
10)  using  a  compound  microscope  with  an  optical  micro¬ 
meter. 


700 

600 

500 

400 

H 

X 

c*  300 

200 

100 

0 


0  20  40  60  80  100  120  140  160  tIO  200 

Time  (minutes) 

Fig.  2.  Schematic  of  turbulent  Reynolds  number  vs.  lime  in 
behavioral  experiments.  Turbulence  became  stationary  in  this  tank 
after  mm  of  grid  oscillation.  Solid  lines  show  the  duration  of 
each  turbulence  treatment,  and  rectangles  indicate  when  the  video 
record  was  taken.  Rest  periods  are  shown  between  turbulence  treat¬ 
ments.  This  example  is  from  replicate  1. 


Flow  velocity  measurements — Flow  velocities  were  mea¬ 
sured  with  a  Dantec  two-axis  laser  Doppler  velocimeter 
(LDV)  after  all  larval  experiments  were  completed,  to  avoid 
any  interference  of  the  laser  or  seeding  particles  with  larval 
behavior.  LDV  is  not  ideal  for  turbulence  measurements  in 
near-zero  mean  flow  (see  “Results”)  but  was  the  best  in¬ 
strument  for  measuring  velocity  remotely.  The  tank  was 
filled  with  filtered  seawater  and  heavily  seeded  with  10-/xm 
hollow  glass  spheres  (density  =  1.05  g  cm"5).  After  10  min 
of  initial  grid  oscillation,  vertical  and  horizontal  velocities 
were  measured  (n  =  2,000-10,000  records)  by  LDV  at  five 
points  in  the  focal  plane  of  the  larval  video  frame  (Fig.  1). 
Flow  measurements  were  made  at  each  turbulence  treatment 
level  in  a  randomized  order  and  then  were  replicated  for  the 
center  point  only. 

Larval  velocity  measurements — Larval  video  records  were 
digitized  using  an  EP1X  PICXJ  SV5  image  capture  board 
with  XCAP-Std  software.  Each  20-min  still-water  video  seg¬ 
ment  was  digitized  into  40  50-frame  sequences  at  2  frames 
s_l  (fps).  The  turbulent  video  segments  were  digitized  into 
30  sequences  per  treatment  replicate  at  5  and  10  fps  for  Rem 
^  200  and  Rem  >  200,  respectively.  The  sequence  length 
was  limited  by  image  buffer  size,  and  capture  intervals  were 
selected  so  that  movements  of  individual  larvae  could  be 
followed  easily  from  frame  to  frame. 

Image  sequences  were  processed  digitally,  to  measure  the 
larval  vertical  velocities.  Larval  centroid  positions  were 
found  using  Matlab  software  (provided  by  S.  M.  Gallager) 
and  linked  into  larval  trajectories  using  a  Matlab  algorithm 
based  on  distance  and  distance/direction  correlations.  Travel 
distances  were  determined  from  a  grid  visible  in  the  video 
frame.  The  grid-spacing  distance  was  calibrated  with  a  ruler 
placed  in  the  focal  plane,  to  account  for  perspective.  1m- 
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probable  trajectory  segments  were  removed  manually.  For 
each  no- flow  replicate,  we  obtained  between  100  and  300 
larval  trajectories.  For  each  turbulence  treatment  replicate, 
we  obtained  several  hundred  to  2,000  trajectories.  To  use 
only  independent  velocity  measurements,  we  took  the  ve¬ 
locity  for  a  single  frame -to- frame  step  from  each  trajectory, 
using  the  step  beginning  closest  to  the  center  of  the  video 
frame.  A  few  outliers  (<2%  of  the  samples)  were  removed 
for  the  statistical  analysis. 

The  software  did  introduce  a  potential  bias  by  splitting 
crossed  trajectories.  Hovering  and  swimming  larvae  move 
more  slowly  than  sinking  larvae,  and  slow-moving  larval 
tracks  had  more  steps  and  a  greater  likelihood  of  being  split. 
The  result  was  an  artificial  increase  in  the  number  of  mea¬ 
surements,  particularly  for  hoverers  and  swimmers.  Because 
this  bias  was  unlikely  to  add  support  to  our  hypothesis,  we 
considered  it  an  acceptable  observation  error. 

Analysis 

Turbulence — To  avoid  disturbance  of  larvae  during  the 
experiments,  we  estimated  turbulence  characteristics  from 
separate  LDV  measurements,  from  relationships  of  grid-tank 
dynamics,  and  from  larval  velocity  measurements,  LDV 
measurements  were  used  to  calculate  the  mean  vertical  flow 
velocity  (w)  and  root-mean-square  (mis)  vertical  velocity  w\ 
using  unbiased  approximations  for  a  burst-sampled  process 
(Buchave  et  al.  1979).  The  spectrum  of  the  velocity  time 
series  was  estimated  with  a  block-averaged,  discrete  perio- 
dogram  for  randomly  sampled  data  (Chan  et  al.  1998).  We 
estimated  the  noise  variance  cr2  and  velocity  variance 
contributions  from  the  spectrum  (Voulgaris  and  Trowbridge 
1998).  The  spectral  rms  velocity  is  Relationships 
of  gnd-iank  dynamics  were  used  to  calculate  the  theoretical 
rms  vertical  velocity  using  the  formulation  in  Hop  finger 
and  Toly  (1976)  with  empirical  constants  given  in  McKenna 
(2000).  This  estimate  of  is  a  function  of  gnd  geometry 
and  forcing  conditions  and  is  linearly  related  to  Reirr.  Larval 
velocity  measurements  were  used  to  estimate  statistically  the 
fluid  velocity  variance  (c2)  during  each  replicate  of  the  lar¬ 
val  experiments.  Maxi  mum- likelihood  estimates  (MLEs)  of 
rms  velocity  (=crj  were  obtained  by  fitting  a  normal 
mixture  model  to  the  larval  velocities,  as  described  in  the 
following  sections. 

For  each  estimate  of  rms  velocity,  the  turbulence  dissi¬ 
pation  rate  e  was  calculated  as  e  =(w')3f -l,  where  the  eddy 
length  scale  £  was  taken  to  be  £  —  0.2z,  and  z  is  the  distance 
from  the  gnd  (Tennekes  and  Lumley  1972;  McKenna  2000). 
The  Kolmogorov  length  scale  was  estimated  by  17  = 
(j/»£- i)02j  where  v  js  thc  kinematic  viscosity  (Tennekes  and 
Lumley  1972). 

Larval  behavior— Our  goal  was  to  determine  how  many 
larvae  were  sinking  with  retracted  velums,  because  sinking 
was  the  only  behavior  that  we  expected  to  be  associated  with 
settlement.  The  proportion  of  sinkers  could  not  be  estimated 
simply  by  integrating  over  the  proportion  of  observed  neg¬ 
ative  (downward)  velocities,  however,  because  the  observed 
larval  velocities  have  both  a  behavioral  component  and  a 
fluid  transport  component,  A  sinking  larva  can  be  carried 


upward  by  an  eddy,  resulting  in  an  observed  velocity  in  an 
upward  (positive)  direction.  Even  in  still  water,  swimming 
and  hovering  larvae  (swimmers  that  maintain  a  relatively 
constant  position)  move  both  upward  and  downward,  mak¬ 
ing  it  impossible  to  classify  larvae  as  sinkers  on  the  basis 
of  downward  movement  alone.  Ideally,  larvae  with  their  vel¬ 
ums  retracted  could  be  classified  visually  as  sinkers,  regard¬ 
less  of  their  velocities.  At  higher  turbulence  levels,  though, 
motion  blur  in  our  video  images  prevented  us  from  classi¬ 
fying  larvae  by  visual  inspection.  Larvae  could  not  be  cat¬ 
egorized  individually  into  behavioral  modes;  instead,  we  an¬ 
alyzed  larval  behavior  statistically,  by  fitting  a  normal 
mixture  model  to  our  larval  velocity  data  to  estimate  the 
proportions  of  swimming,  hovering,  and  sinking  larvae. 

In  our  analysis  of  larval  behavior,  several  assumptions 
were  made. 

(1)  Instantaneous  flow  velocities  Wo  are  normal,  with 
mean  and  variance  c rj. 

W0  —  of) 

(2)  In  still  water,  the  vertical  velocities  W,  of  larvae  in  a 
given  behavioral  mode  arc  approximately  normal,  with 
means  and  variances  a}, 

W,  -  N(jM„crf), 

II  swimmers  (/q  >  0,  21  1) 

2  hoverers  (jq  :S  0,  cr\  ^  1) 

3  sinkers  Qq  <  0,  a]  1,  |jq|  >  jq,  |jq|) 

(3)  Flow  velocity  and  larval  velocity  are  independent  and 
additive — that  is,  relative  larval  velocity  in  flow  WTi  equals 
flow  velocity  W0  plus  larval  behavioral  velocity  W(. 

WT<  =  W„  +  Wt  (1) 

(4)  Larval  sinking  and  swimming  abilities  do  not  change 
significantly  with  flow  conditions.  Within  a  behavioral  mode, 
W ,  is  not  a  function  of  Wa. 

Assumption  1  has  been  validated  in  laboratory  studies  of 
gnd-generated  turbulence  (Mouri  et  al.  2002),  and  our  LDV 
velocity  data  were  normally  distributed.  Assumption  2  was 
justified  by  analysis  of  larval  velocities  in  still  water  (details 
follow).  Assumption  3  is  frequently  used  in  studies  of  par¬ 
ticle  transport  in  turbulence  (e.g.,  Reeks  1977),  because  in¬ 
ertial  forces  are  negligible  for  particles  with  Reynolds  num¬ 
ber  :Sl.  Assumption  4  could  not  be  directly  verified  but  was 
considered  to  be  reasonable  given  the  experimental  condi¬ 
tions.  For  ciliary  swimmers,  the  speed  of  movement  is  di¬ 
rectly  proportional  to  the  ciliary  beat  frequency,  which  is 
limited  by  viscosity  and  by  salinity-  and  temperature-depen¬ 
dent  biochemical  rates  (Chia  et  al.  1984;  Podolsky  and  Emlet 
1993;  Young  1995).  These  seawater  properties  remained 
constant  during  our  experiment,  and  there  would  have  been 
no  viscosity-  or  biochemical  rate-related  effects  on  larval 
swimming  or  hovering  abilities.  Larval  swimming  orienta¬ 
tions  can  potentially  be  affected  by  velocity  shear  (Jonsson 
et  al.  1991),  but  such  an  effect  was  not  apparent  in  our  ex¬ 
periment  (see  “Discussion”  section). 

Our  assumptions  place  no  constraints  on  behavioral 
changes  from  one  mode  to  another.  Within  a  given  behav¬ 
ioral  mode,  the  larval  velocities  are  drawn  from  a  fixed. 
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known  distribution,  but  individual  larvae  can  switch  behav¬ 
ioral  modes  without  restnction. 

Mixture  model  for  larval  behavior  in  still  water — We  first 
estimated  the  velocity  means  p„  variances  of,  and  mixing 
proportions  a,  of  swimmers,  hoverers,  and  sinkers  in  no-ftow 
conditions  by  fitting  a  three -component  normal  mixture 
model  to  the  still-water  larval  velocity  distributions.  The 
probability  density  of  observed  larval  vertical  velocity  wL 
was  modeled  as 


P(wL  |  or„  p„  erf)  =  X  a,P,{wL \p,t  erf) 


(2) 


where,  for  i  =  (1,  2,  3), 


P/K lI  Mr*  o-f) 


-(wL  ~  p,)2 
2af 


(3) 


Maxi  mum- likelihood  parameter  estimates  of  a„  ph  and 
erf  were  calculated  using  the  expectation-maximization  (EM) 
algorithm  (McLachlan  and  Peel  2000). 

In  still  water,  the  proportion  of  sinkers  was  too  low 
(<0.03  in  all  replicates)  to  use  MLEs  of  p^  and  a]  to  de¬ 
scribe  the  sinking  velocities.  The  fall  velocities  of  dead  lar¬ 
vae  were  approximately  normal  and  were  used  as  a  proxy 
for  the  velocity'  distribution  of  Live,  actively  sinking  larvae 
with  retracted  velums.  Although  larvae  can  sink  with  ex¬ 
pended  velums  by  stopping  ciliary  motion,  we  observed  that 
most  larvae  fully  retracted  their  velums  when  sinking  in  tur¬ 
bulence,  Given  the  predominance  of  velum  retraction  and 
the  relatively  small  difference  in  velocities  of  die  two  types 
of  sinkers  (<10%;  H.  L  Fuchs  unpubl.  data),  we  included 
only  sinkers  with  retracted  velums  in  the  model. 

The  variables  p,  and  erf  characterize  the  behavioral  ve¬ 
locities  of  swimming,  hovering,  and  sinking  larvae  from 
each  culture  (Fig.  3a)  and  were  used  as  known  values  (by 
assumption  4)  in  the  analysis  of  larval  behavior  in  turbu¬ 
lence.  Normality  of  the  modes  was  assessed  using  probabil¬ 
ity  plots.  We  also  estimated  the  behavioral  velocity  range  as 
the  difference  between  mean  swimming  and  sinking  veloc¬ 
ities,  calculated  as  AwL  —  pt  —  p^  (Fig.  3a).  Larger  values 
of  Awl  indicate  that  larvae  have  more  control  over  their 
vertical  position  in  the  water. 


Mixture  model  for  larval  behavior  in  turbulence — We  es¬ 
timated  the  proportions  ar<  of  swimmers,  hoverers,  and  sink¬ 
ers  in  each  turbulence  treatment  by  fitting  a  three-compo¬ 
nent,  normal  mixture  model  to  the  observed  velocity 
distributions  of  larvae  in  turbulence.  By  assumption  3,  the 
relative  velocity  of  a  larva  in  turbulence  WTt  is  the  sum  of 
its  behavioral  velocity  W,  and  the  fluid  velocity  W0.  The 
random  variables  WT,  are  distributed  with  normal  probability 
densities  pT>,  with  means  and  variances  given  by 


pr>  ~  M,  +  Me, 
=  <rf  +  a\ 


(4) 


For  each  observed  larval  relative  velocity  wL,  the  behav¬ 
ioral  mode  i  of  the  observed  larva  was  unknown,  and  the 
probability  density  of  wL  in  turbulence  was  modeled  as 


(mm  s  *) 


Fig.  3.  Graphical  representation  of  the  three-component  mixture 
model,  showing  the  contribution  of  swimming,  hovering,  and  sink¬ 
ing  behaviors  to  the  observed  vertical  velocity  distribution  of  larvae 
in  turbulence-  (a)  Behavioral  velocity  distributions  of  sinking,  hov¬ 
ering,  and  swimming  larvae  were  estimated  for  larvae  in  still  water 
(dead  larvae  were  used  to  determine  sinking  velocities)  and  were 
assumed  to  be  independent  of  flow.  Velocity  means  p,  and  SDs  cr, 
are  indicated  for  each  behavioral  mode,  (b)  The  three-component 
mixture  model  for  observed  larval  velocity  distributions  in  turbu¬ 
lence  has  contributions  from  each  behavioral  mode,  where  pr  are 
the  probability  densities  of  individual  modes  and  ar  are  the  mixing 
proportions.  Velocity  means  (/v  =  p„  +  m3  are  indicated  for  each 
behavioral  component,  the  SD  (i/r  =  V(trj  t  erfa)]  is  indicated 
only  for  the  sinking  mode,  and  pa  and  or’  arc  the  fluid  velocity 
mean  and  variance.  When  the  three -component  mixture  model  was 
fitted  to  observed  velocity  data,  ar  and  o0  were  estimated  by  max¬ 
imum  likelihood.  In  this  example,  they  were  chosen  arbitrarily  for 
illustrative  purposes  (arT,  =  0.33,  p„  -  —3.0  mm  s~r,  cr,  =  3.0  mm 
s_l,  and  Rem  ~  100).  Note  that  because  the  mean  flow  velocity  is 
negative,  the  model  predicts  that  77%  of  observed  velocities  would 
be  <0  here,  although  the  proportion  of  sinking  larvae  was  only 
33%. 


P Awl  I  ttrh  Mr'-  <Tr‘)  =  2  <*r'PT'(wjMr'.  &$')  (5) 

(Fig.  3b).  The  mean  flow  velocity  pa  (=(w))  from  LDV  mea¬ 
surements  was  used  as  a  known  value.  By  assumption  4,  p, 
and  erf  from  still-water  segments  were  also  used  as  known 
values  in  eqs.  4  and  5.  The  unknown  parameters  ( aT ,  and 
cr;)  were  estimated  for  each  treatment  replicate  using  a  mod¬ 
ified  EM  algorithm  (Web  Appendix  1  at  http://www.aslo.org/ 
1  o/ toe/ v  oL4  9/ issueJS/ 193 7al.pdf)  to  maximize  the  log-like¬ 
lihood  of  eq.  5  over  all  observed  velocities.  Within  the 
constraints  of  the  model  (eqs.  4,  5),  crr>  has  no  statistical 
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Fig.  4.  Fluid  velocity  characteristics  calculated  from  LDV  measurements,  from  theory  and  from 
observed  larval  velocities,  (a)  Mean  vertical  velocity  (w).  (b)  rms  vertical  velocity  w\  (c)  Turbulence 
dissipation  rate  e.  (d)  Kolmogorov  length  scale  17.  Theoretical  estimates  were  obtained  using  an 
empirical  relationship  between  grid-forcing  conditions  and  rms  velocities  (Hopfinger  and  Toly 
1976).  Direct,  unbiased  estimates  were  calculated  from  raw  LDV  data  (Buchave  ct  al.  1979).  Spec¬ 
tral  estimates  were  calculated  by  removing  the  estimated  noise  contribution  from  LDV  velocity 
spectra  (Voulgaris  and  Trowbridge  1998).  MQUBs  were  obtained  by  fitting  a  three-component  mixture 
mode)  to  observed  velocities  of  larvae  in  turbulence  and  averaging  estimates  from  six  replicates. 
Error  bars  are  1  SE.  Mean  flow  velocity  (a)  cannot  be  estimated  theoretically  or  by  spectral  methods. 


dependence  on  cr*.  and  these  parameters  can  be  estimated 
simultaneously  from  the  observed  larval  velocities 

Relating  behavior  to  turbulence — We  wished  to  determine 
the  general  relationships  between  turbulence  characteristics 
and  the  proportion  of  sinking  larvae.  Maximum- likelihood 
values  of  (  =  £rc  from  the  mixture  model)  were  our  best 
estimates  of  turbulence  intensity  in  individual  replicates  (see 
“Results”  section),  and  these  were  used  to  calculate  the  dis¬ 
sipation  rate  e  and  the  Kolmogorov  length  scale  17  for  each 
turbulence  treatment.  The  parameters  were  averaged  by 
treatment  level,  and  the  relationships  between  e  or  17  and  the 
proportion  of  sinkers  oy 3  were  estimated  by  exponential  re¬ 
gression. 

Any  gTOUp  of  larvae  has  a  net  mean  behavioral  velocity 
(wL),  and  the  relationship  between  (wL)  and  e  is  of  interest 
for  addressing  population-level  questions  about  larval  fluxes. 
We  calculated  the  net  mean  velocity  of  larvae  in  each  treat¬ 
ment  replicate  as  (uy)  =  Z  aTtfjL(.  The  relationship  between 
(wL)  and  log]0£, 

(wy)  -fiy  -  a0e(o,to*l0*)  (6) 

was  estimated  by  exponential  regression.  The  net  velocity 
(w,  )  is  bounded  above  by  the  mean  swimming  velocity  ^ 
and  falls  off  exponentially  with  log10e.  Linear  and  logistic 


functions  of  e  were  also  fitted  to  (wL)  but  had  poorer  fits 
and  are  omitted.  Because  larval  swimming  abilities  varied 
among  cultures,  the  model  was  fitted  for  individual  replicates 
and  for  pooled  data,  excluding  replicates  with  maximum  and 
minimum  values  of  AwL. 

Results 

Turbulence  characteristics — A  comparison  of  LDV  mea¬ 
surements  from  five  points  in  the  video  frame  indicated  that 
turbulence  was  generally  homogeneous  within  our  measure¬ 
ment  area,  although  some  inhomogeueities  appeared  at  the 
highest  turbulence  level.  At  ReiCT  =  600,  w'  values  at  the 
outer  points  differed  from  wr  at  the  center  point  by  5-20%. 
In  spite  of  this  spatial  inhomogeneity,  we  considered  it 
worthwhile  to  present  results  from  all  turbulence  treatment 
levels.  All  reported  LDV  results  refer  to  the  average  of  rep¬ 
licate  measurements  taken  at  the  center  point  of  the  video 
frame. 

The  mean  vertical  flow  velocity  (w)  was  generally  nega¬ 
tive  but  always  had  a  magnitude  <7  mm  s_l  (Fig.  4a).  The 
relationship  between  Rem  and  mean  flow  velocity  appeared 
to  be  nonlinear,  presumably  because  the  pattern  of  circula¬ 
tion  in  the  tank  has  a  nonlinear  dependence  on  forcing  con¬ 
ditions  (McDougaU  1979;  McKenna  2000). 
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Table  3.  Shell  size  and  still-water  velocities  for  competent  larvae,  given  as  mean  i  1  SD.  Fall  velocities  were  measured  for  a  subsample 
of  dead  larvae  (n  =  20)  from  each  replicate.  Velocities  were  measured  for  live  larvae  (100  <  n  <  300)  during  each  replicate,  and  the 
velocities  of  swimmers  and  hoverers  and  the  mixing  proportions  were  estimated  by  maximum  likelihood.  Also  given  is  the  mean  difference 
between  swimming  velocity  and  fall  velocity  (Aw,)  for  each  culture. 


Replicate 

Age 

(days) 

Shell  length 
(^m) 

Fall  velocity 
(mm  s_t) 

Hovering  velocity 
(mm  s  ') 

Swimming  velocity 
(mm  s’1) 

Proportion  of 
sinkers 

Proportion  of 
swimmers 

Awl 

(mm  s_1) 

1 

24 

648±93 

-9.2±  1.1 

-04±07 

4.0+1.2 

003 

0.15 

13.2 

2 

24 

690  ±43 

— 9.2±  1.0 

— 1,3  ±0.2 

1.7±2.9 

001 

0.47 

10.9 

3 

24 

587±40 

— 6.4±  1.2 

-0.1±1.0 

3.0±1.6 

0.00 

0.20 

9.3 

4 

24 

662  ±69 

— 8.8±  1.5 

-0.5  ±0.8 

3.0±1  6 

0.00 

0.25 

11.8 

5 

24 

767±58 

— 9.8±  1.5 

— 0.5±  1.0 

3.6±1.0 

0.00 

0.27 

13.4 

6 

27 

605  ±46 

— 7.0±  1.2 

-0  9±0.3 

0.8±1.8 

000 

0.61 

7.9 

Direct  estimates  of  w'  increased  linearly  with  Rem  ( r2  — 
0.97),  from  10.1  mm  s-1  at  the  lowest  turbulence  level  to 
18.9  mm  s“l  at  the  highest  level  (Fig.  4b).  Spectral  estimates 
of  were  consistently  — 7  mm  s-1  less  than  direct  estimates 
of  w',  indicating  a  substantial  noise  contribution  in  the  mea¬ 
sured  velocity  variance  (Fig.  4b).  The  noise  was  assumed  to 
have  zero  mean  and  no  effect  on  (w).  The  spectral  method 
sometimes  overestimates  the  noise  variance  (Voulgaris  and 
Trowbridge  1998),  but  LDV  measurements  are  generally 
noisy  because  of  Doppler  noise  and  velocity  shear  in  the 
measurement  volume  (Buchave  et  al.  1979;  Voulgaris  and 
Trowbridge  1998).  Because  of  the  high  percentage  of  noise 
in  the  LDV  measurements,  the  direct  and  spectral  calcula¬ 
tions  were  probably  over-  and  underestimates  of  the  rms  ve¬ 
locity,  but  they  are  presented  here  as  independent  reference 
estimates. 

When  it  was  estimated  from  the  three-component  mixture 
model,  the  fluid  velocity  variance  was  always  close  to 
the  theoretical  values  derived  from  empirical  relation¬ 
ships  of  grid  tank  dynamics  (Fig.  4b).  This  was  also  true  for 
the  turbulence  dissipation  rate  e  (Fig.  4c)  and  Kolmogorov 
length  scale  rj  (Fig.  4d),  which  were  calculated  directly  from 
rms  velocity  estimates.  The  close  correspondence  between 
our  average  estimates  of  and  illustrates  the  validity 
of  using  a  mixture  model  to  estimate  turbulence  intensity  by 
maximum  likelihood  from  the  measured  larva]  velocities. 

Because  of  uncertainty  in  w'  and  w'  from  the  LDV  mea¬ 
surements  (due  to  the  noise  contribution)  and  because  MLEs 
of  directly  represent  the  turbulence  intensity  during  each 
replicate  of  the  larval  experiments,  we  discuss  our  larval 
behavior  results  in  terms  of  turbulence  characteristics  cal¬ 
culated  from  MLEs  closely  resembled  theoretical 
values  w'm,  and  the  relationships  between  larval  behavior 
and  are  qualitatively  the  same  as  those  between  behavior 
and  As  calculated  from  w^,  ,  e  was  8.1  X  10“3-2.7  X 
10°  cm2  s-3,  and  -q  was  0.2-L2  mm. 

Larval  velocities  in  still  water — Although  all  larvae  were 
raised  under  identical  conditions,  larval  size  and  swimming 
abilities  varied  between  cultures.  Mean  fall  velocities  were 
6. 4-9. 8  mm  s_I  (Table  3),  and  replicates  3  and  6  had  sig¬ 
nificantly  slower  fall  velocities  than  other  replicates  (one¬ 
way  analysis  of  variance,  Fi  Ui  =  24.268,  p  =  0.001  with 
post-hoc  Tukey’s  test).  The  slowest-sinking  larvae  were  also 


the  smallest,  and  larvae  in  replicate  6  had  the  slowest  mean 
swimming  velocity. 

In  no-flow  conditions,  hovering  and  swimming  larvae 
formed  two  distinct  groups  in  all  replicates  (Fig.  5a).  Swim¬ 
mers  always  had  positive  mean  velocities  and  velocity  var¬ 
iances  >1  mm  s~\  whereas  hoverers  had  negative  mean  ve¬ 
locities  and  velocity  variances  <1.1  mm  $-1  (Table  3).  The 
proportion  of  sinking  larvae  was  always  <0.03  and  was 
<t0.01  in  four  of  six  cultures.  The  proportions  of  swimming 
and  hovering  larvae  varied  between  cultures.  The  behavioral 
velocity  range  A wL  was  7.9-13.4  mm  s"1. 

Larval  behavior  in  turbulence — The  distributions  of  ob¬ 
served  larval  velocities  became  wider  at  higher  turbulence 
intensities  (Fig.  5)  because  of  the  increased  flow  variance. 
The  distributions  also  shifted  toward  more  negative  veloci¬ 
ties  at  higher  turbulence  levels,  in  part  because  of  the  mean 
downward  flow  and  in  part  because  of  the  increasing  pro¬ 
portion  of  sinking  larvae.  For  larvae  in  still  water,  the  ve¬ 
locity  distributions  were  always  positively  skewed  and  bi- 
modal,  whereas  in  turbulence,  the  velocity  distributions  were 
negatively  skewed. 

We  found  a  clear  relationship  between  turbulence  intensity 
and  the  proportion  of  larvae  in  each  behavioral  mode,  as 
estimated  from  the  three -component  normal  mixture  model 
(Fig.  6).  The  proportion  of  sinking  larvae  was  effectively 
zero  at  Rem  —  100  and  generally  increased  with  turbulence 
intensity.  The  proportion  of  hovering  larvae  was  lower  in 
turbulence  than  in  no-flow  conditions.  The  proportion  of 
swimming  larvae  was  higher  than  other  groups  in  turbulence 
up  to  Rem  —  400  (e.g..  Fig.  7)  and  then  decreased  at  higher 
turbulence  levels  as  the  mean  proportion  of  sinkers  in¬ 
creased. 

Relating  behavior  and  turbulence — When  data  were  av¬ 
eraged  by  turbulence  treatment,  the  mean  proportion  of  sink¬ 
ing  larvae  aT,  increased  approximately  exponentially  with 
the  log  of  the  dissipation  rate  [or^  =  0.26  exp(l.l  Iogl0e),  r2 
=  0.89;  Fig.  8a]  and  decreased  approximately  exponentially 
with  the  Kolmogorov  length  scale  [ar  —  0.71  exp(-3.l7}), 
r2  =  0.86;  Fig.  8b]. 

An  exponential  model  provided  a  good  fit  (r1  =  0.73  for 
pooled  data)  to  the  relationship  between  the  net  behavioral 
velocity  (wl)  and  the  dissipation  rate  e  (Fig.  9).  The  net 
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behavioral  velocity  was  negative  at  dissipation  rates  greater 
than  e  ~  7.4  X  10_2-3.2  X  10  1  cm2  s-3.  When  the  data 
were  pooled,  this  threshold  was  at  1.9  x  10_l  era2  s-3. 

Discussion 

Mud  snail  veligers  in  our  study  altered  their  behavior  in 
response  to  turbulent  flow  conditions  and  sank  more  fre¬ 
quently  in  more  intense  turbulence.  Although  we  cannot 
show  conclusively  that  veliger  sinking  behavior  is  a  precur¬ 
sor  to  settlement,  the  response  was  demonstrated  in  late- 
stage,  competent  larvae.  Moreover,  the  turbulence-induced 
response  was  most  pronounced  under  turbulence  conditions 
comparable  to  energetic  nearshore  areas  where  intertidal  spe¬ 
cies  must  settle. 

Our  experimental  turbulence  treatments  were  specifically 
appropriate  for  questions  about  settlement  in  tidal  flows  at 
Barnstable  Harbor,  the  natal  habitat  of  our  larvae.  There  are 
no  published  measurements  of  dissipation  rates  from  the  har¬ 
bor,  but  Hunt  and  Mulljneaux  (2002)  reported  shear  veloci¬ 
ties  of  up  to  =  3.5  cm  s-1  over  the  inner-harbor  tidal 
flats  during  flood  tide.  Using  the  relationship  e  =  m*/kz, 
where  k  is  von  Karman’s  constant  and  z  is  height  above 
bottom  (Gross  and  NoweLl  1985),  the  dissipation  rate  can  be 
estimated  as  £  «  10°  cm2  s-3  over  the  flats.  We  measured 
turbulence  in  Barnstable  Harbor  during  one  complete  tidal 
cycle  in  August  2002  and  estimated  dissipation  rales  to  be 
e  •=»  10-M01  cm2  s'3  (H.  L.  Fuchs  unpubl.  data).  The  lowest 
turbulence  level  used  in  our  experiments  (e  =  8,1  X  10-3 
cm2  s~3)  is  probably  representative  of  slack  tide,  whereas  the 
higher  turbulence  levels  (up  to  e  -  2.7  X  10°  cm2  s'3)  would 
be  found  during  flood  or  ebb  tides.  Although  laboratory- 
generated  turbulence  is  often  more  intense  than  that  in  rel¬ 
evant  field  conditions  (Peters  and  Redondo  1997),  our  ex¬ 
perimental  turbulence  levels  were  comparable  to  those  in  an 
adult  mud  snail  habitat. 

Our  behavioral  results  (Fig.  8)  suggest  that  veligers  will 
exhibit  sinking  behavior  in  turbulent  (e  >  10~2  cm2  s'3), 
shallow  areas  but  not  in  calmer  (e  <  10~2  cm2  s-3)  offshore 
water.  Our  experiments  were  conducted  at  turbulence  levels 
representative  of  tidal  channels  and  partially  mixed  estuaries 
(e  =  I0~3-10°  cm2  s-3;  Table  1),  where  behavioral  changes 
could  potentially  affect  larval  supply  to  intertidal  habitats. 
Dissipation  rates  are  generally  lower  in  offshore  areas  and 
are  unlikely  to  induce  a  larval  sinking  response,  even  during 
moderate  storms.  In  winds  *=15  m  s_l,  dissipation  rales  at 
the  surface  and  thermoctine  reach  £  =  10-3-10~2  cm2  s-3 
(Dillon  and  Caldwell  1980),  but  our  larvae  rarely  sank  at  e 
^  10~2  cm2  s“3.  Although  stormy  conditions  at  sea  could 
cause  some  veligers  to  sink  (Barile  et  al.  1994),  we  expect 
that  larvae  will  encounter  and  respond  to  strong  turbulence 
primarily  in  nearshore  areas. 

Larval  sinking  in  turbulent,  coastal  zones  could  potentially 
affect  horizontal  transport  of  larvae  over  spatial  scales  of 
tens  of  kilometers.  Pringle  and  Franks  (2001)  described  an 
asymmetric  mixing  mechanism  that  causes  sinking  particles 
to  be  transported  shoreward  in  tidal  currents.  It  is  possible 
that  physical  mechanisms  such  as  asymmetric  mixing  trans¬ 
port  could  enhance  the  retention  of  sinking  larvae  in  coastal 
inlets,  providing  additional  opportunities  for  settlement. 
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Fig.  5.  Normalized  histograms  of  observed  larval  vertical  veloc¬ 
ities  for  all  turbulence  levels,  from  replicate  1.  Solid  curves  show 
the  best  fits  from  the  three-component  mixture  model. 


Larval  behavior — We  expected  mud  snail  larvae  to  sink 
or  swim,  but  we  found  that  they  also  hover,  especially  in 
still  water.  Many  species  of  mollusk  larvae  hover  by  pro¬ 
ducing  mucous  strings  that  act  as  natural  tethers  (Fenchel 
and  Ockelmann  2002),  but  our  larvae  had  no  apparent  mu¬ 
cous  strings  and  were  probably  large  and  dense  enough  to 
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Fig.  7.  Normalized  histogram  of  observed  larval  velocities  from 
replicate  1,  Rem  =  200,  showing  fitted  velocity  distribution  for  the 
three-component  mixture  model  (solid  curve)  with  contributions 
from  each  behavioral  mode  (dashed  lines). 
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Fig,  6.  Maxi  mum- Likelihood  values  of  mean  proportions  tty  of 
sinking,  hovering,  and  swimming  larvae  vs.  Re ^  Symbols  show 
averaged  estimates  from  all  replicates  for  each  turbulence  treatment 
level.  Error  bars  are  1  SE. 


hover  without  added  drag.  Bivalve  larvae  create  wider  feed¬ 
ing  currents  and  feed  more  efficiently  by  hovering  than  by 
swimming  (Gallager  1993;  Fenchel  and  Ockelmann  2002). 
In  calm  water,  where  contact  rates  between  larvae  and  food 
particles  are  low,  mud  snail  larvae  also  may  hover  more  to 
increase  feeding  efficiency. 

Larvae  in  our  experiments  sank  more  frequently  at  higher 
turbulence  intensities,  and  this  tendency  could  have  impor¬ 
tant  implications  for  settlement.  Given  constant  turbulence 
conditions,  sinking  particles  reach  the  bottom  more  quickly 
and  accumulate  to  higher  near-bottom  concentrations  than 
neutrally  buoyant  particles  (e  g.,  McNair  et  al.  1997).  Com¬ 
petent  mud  snail  larvae  would  reach  the  benthos  more  quick¬ 
ly  by  sinking  than  by  swimming,  but  there  is  no  apparent 
benefit  for  larvae  that  sink  in  calm,  deep  water  (e  g.,  off¬ 
shore)  over  unsuitable  habitats.  The  exponential  relationship 
between  the  proportion  of  sinkers  and  Iog]0e  suggests  that 
larvae  should  sink  primarily  in  near-shore  areas  during  en- 


Fig.  8.  Mean  proportion  of  sinking  larvae  ar>  vs.  (a)  log10  of  the  dissipation  rate  e  and  (b) 
Kolmogorov  length  scale  77.  All  values  are  from  MLEs.  Symbols  show  average  estimates  from  six 
replicates  per  treatment  level  (excluding  Retn  =  0),  and  error  bars  are  1  $E.  Solid  lines  show 
exponential  regressions. 
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Fig.  9.  Net  mean  behavioral  velocity  of  larvae  (h^)  vs.  logI0  of 
the  dissipation  rate  e.  Symbols  show  velocity  estimates  for  each 
treatment  from  replicates  1-4  (replicates  with  maximum  and  min¬ 
imum  values  of  AwL  were  excluded).  Dashed  lines  show  the  ex¬ 
ponential  regressions  for  individual  replicates,  and  the  solid  line 
shows  the  regression  fit  to  pooled  data,  using  /x,  =  4.1  mm  s"’  and 
Awl  =  13.3  mm  s~'. 


ergetic  tides,  as  might  be  expected  if  the  response  were  re¬ 
lated  to  settlement. 

Veliger  larvae  clearly  respond  to  turbulence,  but  the  phys¬ 
iological  mechanism  of  turbulence  detection  is  unknown. 
Mud  snail  larvae  retract  their  velums  when  the  cilia  are 
touched  (Dickinson  2002),  indicating  the  presence  of  mech- 
anosensory  cilia  that  could  detect  velocity  shear.  It  has  also 
been  hypothesized  that  gastropod  larvae  sense  acceleration 
with  their  statocysts  (Chia  et  al.  1981;  Crisp  1984;  Young 
1995).  Chia  et  al.  (1981)  found  a  neural  connection  between 
the  statocysts  and  cilia  in  nudibranch  veligers  ( Rostanga  pul - 
chrd),  which  suggests  that  these  organs  might  have  related 
functions.  Both  the  velar  cilia  and  the  statocysts  could  po¬ 
tentially  be  used  for  detection  of  shear  and  acceleration  in 
small  eddies. 

It  is  interesting  to  note  that  the  sharpest  increase  in  the 
proportion  of  sinking  larvae  occurred  when  the  Kolmogorov 
length  scale  was  less  than  the  larval  body  length  (Fig.  8b). 
Although  Kolmogorov-scale  eddies  contain  a  very  small  per¬ 
centage  of  the  total  turbulent  kinetic  energy,  velocity  gra¬ 
dients  exist  even  at  scales  smaller  than  millimeters  (Lazier 
and  Mann  1989).  Larvae  may  detect  and  respond  to  shear 
in  the  smallest-scale  eddies,  and  sinking  behavior  could  de¬ 
pend  on  interactions  of  larvae  with  eddies  at  the  Kolmogo¬ 
rov  scale.  Alternatively,  larvae  may  retract  their  velums  in 
response  to  being  accelerated  or  rotated  by  the  flow  vorticity, 
which  increases  with  the  dissipation  rate.  Regardless  of  the 
detection  mechanism,  descriptors  of  small-scale  turbulence 
(p  and  e)  are  probably  the  most  relevant  flow  characteristics 
for  understanding  larval  responses  to  turbulence. 

Because  vorticity  increases  with  dissipation  rate,  velocity 
shear  could  potentially  affect  the  orientation  and  directed- 


swimming  abilities  of  larvae  in  strong  turbulence.  Larvae 
that  normally  swim  in  a  velum-up  orientation  will  begin  to 
tumble  if  the  viscous  torque  caused  by  velocity  shear  across 
the  body  greatly  exceeds  the  gravitational  torque,  caused  by 
an  asymmetric  density  distribution.  Tumbling  is  defined  as 
rotation  past  90°  from  the  normal,  velum- up  orientation 
(Jonsson  et  al.  1991).  We  calculated  the  critical  shear  (as  in 
Kessler  1986;  Jonsson  et  al.  1991)  required  to  tumble  our 
larvae  as  sheara  *=»  102  s_1.  This  is  about  an  order  of  mag¬ 
nitude  greater  than  the  highest  shear  rate  in  our  experiments, 
which  we  estimated  as  shear  15:5  e0-5^-05  =  1.6  X  101  s_1. 
For  shear  <  2  X  10*  s-1,  the  gravitational  torque  on  the 
larvae  will  always  compensate  the  viscous  torque,  allowing 
larvae  to  remain  within  about  10°  of  the  normal  upward 
orientation  (Jonsson  et  al.  1991).  Visual  inspection  of  our 
video  record  confirmed  that,  even  at  the  highest  turbulence 
level,  in-focus  larvae  were  always  oriented  with  the  velums 
up.  We  conclude  that  velocity  shear  had  little  or  no  effect 
on  larval  swimming  abilities  in  our  experiment  and  that 
competent  mud  snail  larvae  probably  rumble  only  in  ex¬ 
tremely  turbulent  conditions  (e  g.,  the  surf  zone). 

Implications  of  turbulence-induced  sinking  behavior — 
There  are  possible  ecological  and  evolutionary  benefits  for 
veligers  that  sink  in  turbulence.  Gastropod  larvae  may  sink 
in  turbulence  as  a  way  of  avoiding  predators  (Young  1995). 
Mud  snail  veligers  responded  to  turbulence  even  in  the  ab¬ 
sence  of  predators,  and  we  infer  that  larval  fluxes  are  af¬ 
fected  by  turbulence-induced  sinking  behavior,  regardless  of 
whether  turbulence  is  generated  by  predators  or  by  physical 
sources.  Abelson  and  Denny  (1997)  suggested  that  hydro¬ 
dynamic  forces  might  alter  larval  behavior  and  even  provide 
a  settlement  cue.  We  hypothesize  that  a  turbulence- induced 
sinking  response  enables  larvae  to  move  toward  the  bottom 
when  they  reach  shallower  coastal  waters.  Near-bottom  lar¬ 
vae  contact  the  bottom  more  frequently  than  those  in  the 
upper  water  column  (McNair  et  al.  1997)  and  could  have 
more  opportunities  to  test  substrates  and  settle  in  suitable 
areas. 

To  understand  biophysical  coupling  between  turbulence 
and  settlement,  it  is  necessary  to  determine  the  hydrodynam¬ 
ic  conditions  where  larval  behavior  can  affect  sinking  fluxes 
(Crimaldi  et  al.  2002).  The  net  behavioral  velocity  of  our 
larvae  shifted  from  positive  (upward)  to  negative  (down¬ 
ward)  when  the  dissipation  rate  reached  *=  10"1  cm2  s-5. 
This  shift  could  result  in  significant  changes  in  the  advective 
component  of  larval  vertical  fluxes.  On  the  other  hand,  the 
observed  velocity  distributions  of  our  larvae  grew  wider  as 
w'  increased,  indicating  that  larval  movement  is  dominated 
more  by  flow  as  turbulence  intensifies.  Larvae  probably  de¬ 
tect  eddies  at  or  near  the  Kolmogorov  scale,  which  is  a  func¬ 
tion  of  e,  and  we  expect  the  proportion  of  sinking  larvae 
aT>  to  increase  exponentially  with  log10e  in  the  field  as  in 
the  lab.  Behavioral  changes  may  affect  sinking  fluxes  less, 
however,  when  the  vertical  turbulence  intensity  w'  greatly 
exceeds  the  behavioral  velocity  range  AwL, 

Although  we  focused  on  water-column  processes,  it  is 
worth  noting  that  settlement  fluxes  depend  on  both  the  prob¬ 
ability  of  reaching  the  bottom  and  the  probability  of  attach¬ 
ment  to  substrates  (Gross  et  al.  1992;  Crimaldi  et  al.  2002). 
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Like  the  dissipation  rate,  Reynolds  stress  increases  with 
shear  velocity,  and  those  flow  conditions  that  induce  larval 
sinking  are  also  more  likely  to  cause  bed-load  transport  or 
saltation  of  sediments  and  to  prevent  larvae  from  attaching 
to  the  bottom.  In  energetic  flows  such  as  the  strong  tides  at 
Barnstable  Harbor,  a  large  proportion  of  competent  larvae 
could  be  sinking  at  a  given  time  (Fig.  8),  concentrating  near¬ 
er  to  the  bottom  than  nonsinking  larvae  do.  Although  the 
probability  of  larval  attachment  to  substrates  is  lower  in  the 
higher  shear  stresses  of  turbulent  flows  (Gross  et  al.  1992; 
Pawlik  and  Butman  1993),  near- bottom  larvae  could  have 
more  opportunities  to  test  substrates  and  burrow  into  sedi¬ 
ments  during  slack  tides  (tens  of  minutes)  or  brief  lulls  in 
turbulence  (seconds  to  minutes;  Crimaldi  et  al.  2002). 

Turbulence  varies  spatially  and  temporally  in  coastal  ar¬ 
eas,  and  large-scale  settlement  patterns  could  be  influenced 
by  flow -mediated,  active  settlement  processes  rather  than  by 
passive  deposition  alone.  Our  motivation  is  to  understand 
how  biophysical  coupling  between  turbulence  and  settlement 
behavior  affects  the  supply  of  gastropod  larvae  to  coastal 
populations.  We  have  shown  that  turbulence  alters  the  be¬ 
havior  of  competent  mud  snail  larvae  in  the  laboratory.  Field 
and  modeling  studies  are  in  progress  to  determine  how  this 
biophysical  coupling  affects  larval  supply  and  settlement 
fluxes  of  intertidal  gastropods. 
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Chapter  3 


Sinking  in  turbulence:  Effects  of 
larval  behavior  on  larval  supply 
and  settlement  in  tidal  currents 

3.1  Introduction 

Gastropod  larvae  retract  their  velums  and  sink  in  strong  turbulence,  and  it  is  reasonable  to 
hypothesize  that  this  response  is  selected  for  because  it  confers  some  advantage  to  individual 
larvae.  Mud  snail  larvae  ( Ilyanassa  obsoleta )  have  three  behaviors  (swimming,  hovering, 
and  sinking),  and  the  proportion  of  larvae  engaged  in  these  behaviors  depends  on  the 
turbulence  dissipation  rate  e  [10].  The  dissipation  rate  is  the  rate  of  transfer  of  turbulent 
kinetic  energy  to  heat  energy  at  the  smallest  eddy  scales,  and  is  a  relevant  statistic  for 
quantifying  turbulence  at  the  larval  scale.  Mud  snail  larvae  sink  more  frequently  with 
increasing  dissipation  rate,  resulting  in  a  clear  shift  in  the  average  larval  velocity  from 
upwards  to  downwards  at  e  «  10“]  cm2  s“3  [10].  The  question  is  how  do  individual  larvae 
benefit  from  this  behavior? 

Two  hypotheses  have  been  put  forth  to  explain  why  gastropod  larvae  sink  more  fre¬ 
quently  in  turbulence.  The  first  is  that  larvae  of  coastal  gastropods  use  turbulence  as  an 
initial  settlement  cue  [3,  10].  I  use  “cue”  to  mean  an  indicator  of  potentially  favorable  habi¬ 
tats1.  Turbulence  dissipation  rates  of  >  lO"*1  cm2  s-3  are  uncommon  in  shelf  regions  or 
open  ocean  [11,  6,  23],  but  are  typical  of  coastal  areas  and  tidal  inlets  [12,  14,  29].  Sinking 


lltcue.”  Def.  2.  Merriam  Webster’s  Collegiate  Dictionary.  11th  ed.  2003. 
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in  strong  nearshore  turbulence  would  allow  larvae  to  reach  bottom  more  quickly  and  explore 
settlement  sites  in  potentially  favorable  habitats.  The  second  hypothesis  is  that  gastropod 
larvae  sink  in  turbulence  to  avoid  predators  [31].  Swimming  predators  of  all  sizes,  including 
krill,  herring,  tuna,  and  killer  whales,  generate  turbulence  with  dissipation  rates  estimated 
on  the  order  of  lO-1  cm2  s-3  [19].  Sinking  in  turbulence  at  or  above  this  level  could  pro¬ 
vide  an  escape  mechanism  for  larvae  that  have  no  other  defense  against  predators.  Note 
that  predator-generated  turbulence,  nearshore  turbulence,  and  the  turbulence  threshold  for 
larval  sinking  all  share  an  approximate  lower  limit  of  a  «  10-1  cm2  s-3.  Given  this  com¬ 
mon  turbulence  threshold,  neither  hypothesis  can  be  rejected  based  on  current  empirical 
evidence. 

I  use  a  model  to  explore  the  settlement-cue  hypothesis,  by  evaluating  the  theoretical 
effects  of  turbulence-induced  sinking  on  larval  supply  and  settlement  success  in  tidal  cur¬ 
rents.  In  this  paper  I  define  larval  supply  as  the  concentration  of  larvae  within  1  cm  of 
the  bed,  settlement  success  as  the  probability  that  a  larva  will  settle  within  one  tidal  cycle, 
and  settlement  as  permanent  attachment  to  the  bottom.  Turbulence-induced  sinking  must 
affect  both  the  spatial  and  temporal  patterns  of  larval  supply  and  settlement,  but  this  pa¬ 
per  addresses  only  temporal  patterns  that  are  due  to  tidal  variation  in  turbulence.  Larvae 
that  sink  rapidly  in  turbulence  are  expected  to  have  greater  sinking  fluxes  than  larvae  with 
near-neutral  buoyancy,  particularly  during  flood  and  ebb  tides.  Therefore  the  supply  of 
larvae  to  the  bed  is  expected  to  show  different  periodicity  for  behaving  vs.  passive  larvae. 

The  timing  of  larval  contact  with  the  bottom  could  have  significant  effects  on  settlement 
success,  particularly  in  strong  tidal  currents  where  slack  tides  are  of  short  duration.  In  order 
to  settle,  larvae  must  attach  or  burrow  into  the  substrate,  and  their  ability  to  do  so  is  affected 
by  time-dependent  near-bed  hydrodynamics  [5].  During  flood  and  ebb  tides  when  bed  shear 
stress  is  above  some  critical  value,  sediment  and  larvae  are  transported  as  bedload.  When 
both  the  bed  and  the  larvae  are  mobile,  larval  attachment  to  the  bottom  is  less  probable. 
During  slack  tides  larvae  may  settle  more  easily,  but  the  period  when  shear  stress  is  below 
a  critical  value  can  be  as  short  as  10  minutes  [1].  Temporal  patterns  of  both  larval  supply 
and  attachment  probability  must  be  considered  in  evaluating  whether  turbulencc-induced 
sinking  could  be  a  beneficial  settlement  strategy. 

The  first  goal  of  this  chapter  was  to  develop  predictions  about  the  supply  and  set¬ 
tlement  of  gastropod  larvae  with  turbulence-induced  sinking  behavior  in  tidal  inlets.  A 
1-dimensional,  vertical  advection-diffusion  (AD)  model  was  used  to  characterize  these  dy¬ 
namics  for  passive  and  behaving  larvae.  In  this  model  the  vertical  diffusion  of  larvae  is 
due  to  turbulence  as  specified  through  a  simple  eddy  diffusivity  function,  and  vertical  ad- 
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vection  is  due  to  larval  behavior  as  specified  by  a  turbulence-dependent  vertical  velocity. 
The  AD  model  has  been  used  to  describe  settlement  dynamics  of  passive  larvae  in  steady 
currents  over  variable  roughness  elements  [9],  and  of  passive  [15]  and  negatively  phototactic 
[8]  larvae  in  tidal  currents.  Although  complex  behaviors  can  be  modeled  only  implicitly  by 
AD,  this  remains  the  most  useful  spatially-explicit  model  for  the  movement  of  organisms  in 
environmental  gradients  [16]. 

The  second  goal  of  this  chapter  is  to  determine  how  accurately  the  transport  and  set¬ 
tlement  dynamics  of  larvae  with  multiple  individual  behaviors  can  be  captured  by  an  AD 
model  with  implicit  behavior.  Mud  snail  larvae  have  three  distinct  behaviors,  each  with 
a  different  velocity  distribution.  The  shapes  of  these  distributions  could  affect  the  overall 
settlement  dynamics  in  ways  that  cannot  be  captured  by  a  mean-field  AD  approach.  To 
characterize  any  shortcomings  of  the  AD  model,  an  equivalent  particle-tracking  (PT)  model 
[e.g.,  7,  25]  was  used  to  simulate  the  movement  of  larvae  with  explicit  individual-level  be¬ 
haviors.  Individual-based  PT  models  allow  greater  flexibility  in  defining  rules  of  behavior, 
but  require  more  computation  time  than  concentration-based  AD  models. 


3.2  Model  Descriptions 

The  1-dimensional  vertical  advection-diffusion  model  and  particle- tracking  model  are  out¬ 
lined  with  details  about  solution  methods  in  Sections  3.2.1  and  3.2.5.  The  parameters  that 
define  the  hydrodynamics,  larval  behavior,  and  settlement  in  the  models  are  described  in 
more  depth  in  Sections  3.2.2,  3.2.3,  and  3.2.4. 


3.2.1  Advection-Diffusion  Model 


I  modeled  the  vertical  movement  and  settlement  of  larvae  in  a  turbulent  tidal  boundary 
layer  using  the  1-dimensional  vertical  advection-diffusion  model, 


OC 

dt 


d_ 

dz 


„  „dC 

wC-K  — 


(3.1) 


The  solution  of  (3.1)  predicts  the  concentration  of  larvae  C(z,t)  at  height  z  and  time  t. 
In  this  model  larvae  are  advected  by  their  behavioral  velocity  w(zyt)  and  diffused  by  the 
turbulent  eddy  diffusivity  Larvae  were  modeled  as  passive  (constant  w)  or  with 

turbulence-dependent  behavior,  where  w  varies  with  the  turbulence  dissipation  rate  e{z,t) 
(Section  3.2.3).  Both  the  dissipation  rate  and  the  eddy  diffusivity  K  are  functions  of  depth 
and  a  tid ally-oscillating  shear  velocity  u*,  which  is  proportional  to  the  free-stream  velocity 
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Uh  (Section  3.2.2). 

The  boundary  conditions  of  (3.1)  are 


—wC  +  K  —  0  at  z  =  H 

—wC  +  =  —  $s  at  z  =  0. 


(3.2) 


There  is  no  flux  at  the  surface  ( z  =  i7)  because  larvae  cannot  leave  through  the  air-water 
interface.  A  settlement  flux  is  specified  at  the  bottom  (z  =  0)  as 


$3W  =  ^)C(o,t) 


(3.3) 


where  s  is  a  settlement  velocity.  In  this  model  the  potential  for  substrate-selection  behav¬ 
ior  is  ignored,  and  larvae  are  allowed  to  settle  only  when  they  contact  the  bottom.  The 
settlement  velocity  was  modeled  as  zero,  constant,  or  a  function  of  shear  velocity  (Section 
3.2.4). 

With  time-  and  space- varying  advection  and  diffusion  coefficients,  (3.1)  had  to  be  solved 
numerically.  The  partial  differential  equation  was  solved  with  the  Matlab  6.5  PDE  solver, 
which  discretizes  the  spatial  components  of  the  equation  to  generate  an  ODE  in  time  that 
is  solved  using  numerical  integration.  The  equation  was  solved  on  a  linearly-spaced  grid  of 
1-cm  depth  increments  from  0  to  H  for  H  =  [100, 300,  500, 700]  cm,  using  maximum  current 
speeds  of  =  [10, 35, 65]  cm  s-1.  Time  intervals  of  100  s  were  used  over  one  semidiurnal 
tidal  cycle  (12.25  h),  with  a  uniform  initial  distribution  and  a  —1  hour  spin-up  time.  Tidal 
mixing  quickly  overcame  the  initial  distributions,  and  no  settlement  was  allowed  during  the 
spin-up  time. 


3.2.2  Physical  Parameters 

Mud  snails  live  primarily  in  soft-substrate  intertidal  areas,  and  this  model  was  intended  to 
simulate  larval  settlement  in  a  shallow  tidal  channel.  I  use  Barnstable  Harbor,  MA  (USA) 
as  a  reference  tidal  inlet  because  mud  snails  are  abundant  at  this  site.  The  model  is  run  for 
water  depths  ( H  =  [100,  300,  500,  700]  cm)  that  are  representative  of  the  harbor  at  various 
locations  and  tidal  stages.  Three  maximum  current  velocities  (t/#^  =  [15,35,65]  cm  s"”1) 
were  used,  the  fastest  of  which  is  representative  of  Barnstable  Harbor  [1,  Ch.  4];  the  two 
slower  current  velocities  can  be  taken  to  represent  calmer,  more  sheltered  estuaries.  The 
boundary  layer  is  assumed  to  be  depth-limited,  as  is  the  case  at  Barnstable  Harbor  under 
most  conditions  [1,  Ch.  4]. 
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Tides  are  modeled  as  symmetric,  and  the  along-current  free-stream  velocity  U n  varies 
periodically  as 

U„{t)  =  0.5^^1-008^,  (3.4) 

where  Uh ^  is  the  maximum  free  stream  velocity,  H  is  the  water  depth,  and  T  is  a  tidal 
period  of  12.25  h.  The  shear  velocity  n*  influences  larval  settlement  velocity  (section  3.2.4), 
and  is  related  to  the  free-stream  velocity  by 


u+{t) 


Unit) 

10 


(3.5) 


This  estimate  is  based  on  the  drag  coefficient  relationship  Cd  =  ul/Ufj  [14].  Cj  is  typically 
on  the  order  of  3  x  10-3  for  flow  over  smooth  sandy  substrates  [13,  17],  but  was  found  to  be 
Cd  ~  10-2  in  Barnstable  Harbor,  probably  due  to  form  drag  over  large  sandwaves  (Chapter 

4)- 


The  vertical  eddy  diffusivity  K  controls  turbulent  mixing  in  the  model,  and  is  related 
to  shear  velocity  by 

K  =  u+kz  (3.6) 

where  k  —  0.4  is  von  Karman’s  constant  and  £  is  height  above  the  bottom  [14]  (Fig.  3-la). 
This  linear  form  of  eddy  diffusivity  assumes  that  the  mixing  length  is  proportional  to  depth 
[e.g.,  21],  and  is  valid  near  the  bottom  (i.e.  in  the  log  layer)  but  less  accurate  than  other 
turbulence  closure  schemes  in  the  upper  water  column.  Nevertheless,  (3.6)  should  provide 
a  good  estimate  of  larval  concentrations  in  the  lower  water  column  [15].  A  minimum  value 
of  K  ==  1  was  used  to  prevent  numerical  instabilities  in  the  advection-diffusion  model. 


The  diffusivity  term  parameterizes  large-scale  turbulence  in  the  model,  but  larvae  must 
detect  and  respond  to  small-scale  turbulence.  The  smallest  eddies  are  characterized  by  the 
Kolmogorov  length,  time,  and  velocity  scales,  all  of  which  can  be  expressed  as  functions  of 
the  turbulence  dissipation  rate  e  [28].  In  this  model  s  determines  larval  behavior  (Section 
3.2.3),  and  is  related  to  shear  velocity  by 


€  = 


(3.7) 


(Fig.  3-1  b). 
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U  =  15  cm  s  1  U„  =  65  cm  s 

Hoo  Woo 


K  (cm2  s  X)  K  (cm2  s"1) 


0l_l - * -  ii"  ' 

-1  -0.5  0  0.5 

w  (cm  s"1) 


logj0B  (cm2  s  3) 


w  (cm  s  l) 


Figure  3-1:  Examples  of  modeled  (a-b)  eddy  diffusivity  if,  (c-d)  turbulence  dissipation 
rate  e1  and  (e-f)  larval  velocity  w  =  f{e)  vs.  height  above  bottom  z.  Lines  show  15-min 
intervals  from  slack  tide  to  flood/ebb  tide  for  H  —  5  m  and  two  flow  speeds  U - 
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3.2.3  Behavioral  Parameters 


Larvae  were  modeled  as  passive  or  as  changing  their  behavior  in  response  to  turbulence. 
Passive  larvae  were  modeled  with  a  constant  velocity  of  w  =  0.05  cm  s_1  or  w  ==  —0.05 
cm  s~\  representing  positive  and  negative  buoyancy  respectively.  These  values  are  on  the 
order  of  the  vertical  velocities  reported  for  some  bivalve  and  polychaete  larvae  [e.g.,  4,  20,  2]. 
These  constant  velocities  can  also  be  considered  as  subsets  of  mud  snail  larval  behaviors; 
w  =  —0.05  cm  s-1  is  about  the  average  velocity  of  hovering  larvae,  and  w  =  0.05  cm  s_1  is 
about  the  average  velocity  of  hovering  and  swimming  larvae  combined  in  still  water  [10]. 

Behaving  larvae  were  modeled  as  having  3  behavioral  modes:  1)  swimming,  2)  hovering, 
and  3)  sinking  (e.g.,  Ch.  2  Fig.  3  a).  The  proportion  ct{  of  larvae  engaged  in  mode  i 
depends  on  the  turbulence  dissipation  rate,  as  determined  by  fitting  the  following  functions 
to  laboratory  data  (Fig.  3-2  a-c): 


ai 

=  h¥) 

=  1  —  ai  —  c*3 

(3.8) 

C*3 

=  /3(e) 

where 

'£•)  =  - • 

1  +  exp(-6i0  -  bi i  log10  e) ' 

The  population-average  vertical  velocity  of  larvae  at  a  point  z  in  time  t  is 

3 

w(ztt)  = 

1=1 

Here  /i*  is  the  mean  vertical  velocity  of  larvae  in  mode  i  (Table  3.1).  The  larval  velocity  w 
decreases  as  the  dissipation  rate  increases  (Fig.  3-2  d).  As  a  result,  larvae  on  average  sink 
more  near  the  bottom  than  near  the  surface,  and  sink  more  during  flood  and  ebb  tides  than 
during  slack  tides  (Fig.  3-1  e,  f). 


(3.9) 


(3.10) 


3.2.4  Estimating  Larval  Supply  and  Settlement 

To  characterize  the  effects  of  behavior  on  temporal  patterns  of  larval  supply,  the  AD  model 
was  run  with  no  settlement.  Larval  supply  was  calculated  as  the  concentration  in  the  bottom 
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a 


logiQe(cm2  s  3) 

Figure  3-2:  Proportions  a*  of  (a)  swimming,  (b)  hovering,  and  (c)  sinking  larvae  vs. 
dissipation  rate  logio£  and  (d)  average  larval  vertical  velocity  w  vs.  dissipation  rate  £. 
Circles  are  estimates  from  laboratory  experiments  [10],  solid  lines  are  the  form  used  in 
model. 
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Table  3.1:  Parameters  used  in  the  turbulence-dependent  behavioral  velocity  functions.  The 
vertical  velocity  for  larvae  with  behavior  i  is  shown  as  mean  ±  1  SD.  The  proportions  a,  of 
larvae  engaged  in  each  behavior  are  functions  of  the  dissipation  rate,  shown  in  Fig.  3-2  a-c. 


Behavior  i 

Velocity  (cm  s 

Hi  ±  CT; 

(1)  swimming 

0.41  ±  0.16 

(2)  hovering 

-0.05  ±  0.08 

(3)  sinking 

-0.92  ±  0.12 

1  cm  normalized  by  the  total  number  of  larvae, 


Cb(t)/CT(t)  = 


Jo  C(z,t)dz 


(3.11) 


To  explore  how  the  timing  of  larval  supply  affects  settlement  success,  larvae  in  contact 
with  the  bottom  were  allowed  to  attach  and  settle  with  a  settlement  velocity  s.  Once 
settled,  larvae  were  unable  to  re-enter  the  water  column.  Three  settlement  functions  were 
used  (Fig.  3-3): 

s(t)  =  C\  (constant) 


s(t)  —  C2  (l  U* {£) / IL+max) 

\  0.01  if  u*(t)  <  u*c 
s(t)  =  < 

I  0  if  u*(t)  >  u*c 


(linear) 

(step) 


(3.12) 


The  Linear  and  step  cases  include  a  dependence  on  the  tidally-oscillating  shear  velocity, 
where  u*ma:c  =  0.1  and  u ^  is  the  critical  shear  velocity  for  bedload  transport.  The 
linear  and  step  cases  are  more  realistic  than  a  constant  settlement  velocity,  because  larvae 
are  expected  to  have  more  difficulty  attaching  to  the  bottom  at  higher  shear  velocities.  The 
u+cr  for  mud  snail  larvae  was  assumed  to  be  equal  to  that  of  average  Barnstable  Harbor 
sediment  (diameter  =  100  /zm  [26,  27]),  estimated  from  a  Shields  diagram  to  be  u*cr  ~  1.2 
cm  s-1. 


In  the  step  function,  larvae  have  a  small  settlement  velocity  ( s  =  0.01)  when  the  shear 
velocity  is  below  the  critical  value.  In  the  constant  and  linear  functions,  C\  and  o 2  were 
selected  so  that  the  settlement  velocity  averaged  over  a  tidal  cycle  was  the  same  as  for  the 
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Figure  3-3:  Settlement  velocities  s  (3.12)  vs.  time  for  three  values  of  U .  Line  colors 
indicate  settlement  function:  constant  s ,  light  gray;  linear  s,  medium  gray;  step  s ,  black. 


step  function, 


s  = 


1  fT 

T  L  s^t')dL 


(3.13) 


The  time-averaged  settlement  velocity  was  greater  for  model  runs  with  a  smaller  maximum 
current  velocity  (s  —  0.0066  at  Ujj^  =  15;  s  =  0.0038  at  U =  35;  $  —  0.0027  at  U //^  = 
65),  because  u*(t)  <  u+cr  for  longer  time  periods.  It  is  reasonable  for  the  settlement  velocity 
to  be  higher  in  slower  flows,  because  slower  flows  are  less  likely  to  transport  sediments  and 
exert  weaker  drag  forces  on  larvae  as  they  are  trying  to  settle. 

Settlement  success  Q  was  characterized  by  the  proportion  of  larvae  that  settled  within 
one  tidal  cycle, 

Jo  C(z,T)dz 


Q  =  1 


f*Ci(z)dx 


(3.14) 


where  Cv  is  the  initial  concentration  distribution. 


3.2.5  Particle- Tracking  Model 

In  the  AD  model,  multiple  larval  behaviors  are  implicit  in  the  dissipation-dependent  larval 
velocity  w  (Section  3.2.3).  To  determine  if  an  implicit  behavior  model  accurately  cap¬ 
tures  the  dynamics  of  larvae  with  multiple  behaviors,  larvae  were  also  modeled  using  an 
individual-based  particle-tracking  (PT)  approach.  This  model  allowed  behaviors  to  be  mod- 
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eled  both  implicitly,  as  in  the  AD  model,  and  explicitly  with  stochastic  behaviors  for  each 
individual.  The  position  zn(t  -b  At)  of  larva  n  at  time  t  +  At  is  given  by 


zn(t  +  At)  =  zn(t)  + 


Wn  + 


dK_ 

dz 


A  t  + 


_ _  i  .  ,  IdK A  i 

2K  [zn(t)  +  ~— At,  t)  At 


0.5 


(3.15) 


deterministic 


random 


This  is  the  random-walk  equivalent  of  the  1-D  advection  diffusion  equation,  derived  from 
the  Fokker-Planck  equation  [7].  The  PT  model  is  corrected  for  spatially- varying  diffusivity 
[18,  30,  25],  but  this  correction  is  unnecessary  with  linear  diffusivity.  The  second  term  on 
the  right  hand  side  of  the  equation  gives  the  larva  a  deterministic  movement,  depending  on 
the  larval  velocity  wn  and  the  gradient  in  eddy  diffusivity  dK/dz  at  position  zn(t)  and  time 
t.  The  last  term  of  the  equation  gives  the  larva  an  additional  random-diffusive  movement, 
where  &  ~  Ar(0,l).  In  (3.15)  larval  behavior  is  implicit  and  deterministic,  as  in  the  AD 
model.  When  larval  behavior  is  modeled  explicitly,  the  behavior  term  is  stochastic  (Section 
3.2.3). 

Because  of  the  random  term  in  (3.15),  larvae  were  sometimes  transported  out  of  the 
water  column  to  zn(t  +  At)  >  H  or  <  0.  No-flux  boundary  conditions  were  enforced  by 
reflecting  the  larvae  back  using 


zn{t  At) 


—zn(t  -h  At)  if  zn(t  -f  At)  <  0 

2 H  —  zn(t  +  At)  if  zn(t  4-  At)  >  H. 


(3.16) 


For  comparison  against  the  AD  model,  the  PT  model  was  run  with  H  —  100  for  passive 
larvae  with  three  different  velocities  (wn  —  —0.05,  -0.2,  and  -0.5  cm  s'-1)  and  for  behaving 
larvae  with  deterministic  and  stochastic  behavior  (Section  3.2.3),  always  with  no  settlement. 
The  advection  diffusion  model  allowed  behaviors  to  be  represented  only  as  a  population- 
averaged,  deterministic  velocity  (3.10),  whereas  the  particle  tracking  model  enabled  the  use 
of  individual  stochastic  behaviors  for  each  larva.  The  PT  model  was  run  with  two  levels 
of  complexity.  First,  larval  velocity  was  calculated  as  a  population  average  (3.10)  as  in  the 
AD  model.  Second,  each  larva  was  assigned  to  behavior  i  with  probability  a*  at  each  time 
step,  and  the  velocity  of  each  larva  in  behavior  i  was  drawn  from  a  normal  distribution  with 
mean  /ix  and  variance  (Table  3.1), 


P  {wn(zn(t),t)  ~  =  ai(zn(t),t). 


(3.17) 


The  PT  model  was  run  for  a  complete  tidal  cycle  with  a  ~l-hour  spin-up  time  and  a 
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uniform  initial  distribution.  Four  different  time  steps  (At  =  0.125,  0.25,  0.5,  and  1.0  s)  were 
used  to  characterize  the  errors  associated  with  time-step  size.  PT  solutions  were  taken  as 
the  average  of  10  model  runs.  Because  the  time  steps  differed  in  the  AD  and  PT  models, 
PT  solutions  were  averaged  over  100  s  intervals  so  that  solution  time  points  corresponded 
to  those  of  the  AD  model.  For  each  model  comparison,  model  outputs  were  judged  by  the 
similarity  of  larval  concentration  over  all  z  and  t  grid  points. 

The  AD  and  PT  models  include  two  main  simplifying  assumptions  about  behavior.  First, 
responses  of  larvae  to  turbulence  are  the  only  behaviors  included  in  the  model,  even  though 
real  mud  snail  larvae  exhibit  sediment-selection  behavior  [27]  that  probably  influences  the 
distribution  of  settled  larvae.  In  the  model,  settlement  requires  contact  with  the  bottom  - 
there  is  no  “reaction  distance”  [e.g.  15]  over  which  larvae  can  respond  to  the  bottom,  because 
there  is  no  settlement-related  behavior.  Larval  “tumbling”  [e.g.  20,  24]  is  also  omitted  - 
larvae  that  fail  to  settle  on  contact  with  the  bed  are  not  allowed  to  remain  on  the  bottom  by 
rolling.  The  second  assumption  is  that  behaviors  depend  only  on  instantaneous  dissipation 
rates  and  not  on  the  larval  history  of  behavior.  For  example,  in  the  PT  model  larvae  can 
switch  behavioral  modes  from  one  time-step  to  the  next,  even  though  real  larvae  probably 
have  some  minimum  response  time  before  they  can  switch  behaviors.  These  behavioral 
simplifications  make  the  models  more  tractable  and  make  it  easier  to  evaluate  the  effects  of 
turbulence-induced  sinking  on  larval  supply  and  settlement. 
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Figure  3-4:  Larval  supply  vs.  time  for  three  values  of  £7//^.  Near-bed  larval  concentrations 
Cb  are  normalized  by  the  depth-integrated  concentration  Cy.  Line  colors  indicate  behavior: 
w  —  0.05  cm  s”1,  light  gray;  w  =  —0.05  cm  s_1,  medium  gray;  w  —  /(e),  black.  H  =  5  m. 

3.3  Results 

3.3.1  Larval  Supply 

Temporal  patterns  of  larval  supply  are  very  different  for  passive  larvae  than  for  those  with 
turbulence-dependent  behavior  (Fig.  3-4).  For  all  current  speeds,  negatively-buoyant  larvae 
(w  =  —0.05  cm  s-1)  have  large  peaks  in  near-bed  concentrations  at  slack  tides.  These  peaks 
grow  slightly  smaller  and  narrower  as  the  maximum  current  velocity  increases.  In  contrast, 
positively  buoyant  larvae  (re  —  0.05  cm  s"1)  have  elevated  near-bed  concentrations  during 
flood  and  ebb  tides.  These  “peaks”  are  actually  very  low  near-bed  concentrations  however, 
and  they  change  little  with  increases  in  the  maximum  current  velocity.  The  supply  of 
behaving  larvae  (w  =  /(e))  to  the  bed  is  also  greatest  during  flood  and  ebb  tides,  but  is 
strongly  dependent  on  the  maximum  current  velocity.  For  U //^  =  15  cm  s-1,  the  dissipation 
rate  is  always  low  enough  that  the  larvae  are  swimming  upward  most  of  the  time  (Fig.  3-le); 
very  few  of  these  larvae  concentrate  near  the  bed.  For  =  35  cm  s-1,  behaving  larvae 
have  broad  peaks  in  near-bed  concentration  during  flood/ebb  tides.  For  U =  65  cm  s_1 
these  peaks  are  even  broader,  but  there  is  a  slight  reduction  in  near-bed  concentration  at 
peak  flood  and  ebb  tides  because  of  intense  turbulent  mixing.  For  —  35  and  65  cm 
s'1  ,  behaving  larvae  have  higher  near-bed  concentrations  than  negatively  buoyant  larvae 
during  most  of  the  tidal  cycle,  excepting  slack  tides. 
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Figure  3-5:  Proportion  of  larvae  settled  in  1  tidal  cycle  Q  vs.  water  depth  H  for  three 
flow  velocities  UHoo  (a,  d,  g,  UHoo  =  15  cm  s-1;  b,  e,  h,  UHoo  =  35  cm  s-1;  c,  f,  i,  UJJoo  =  65 
cm  s-1)  and  three  settlement  velocity  functions  s  (a,  b,  c,  constant  s;  d,  e,  f,  linear  s;  g,  h, 
i,  step  s).  Symbols  indicate  behavior  function:  w  —  0.05  cm  s-1,  A;  w  —  —0.05  cm  s_1,  V? 
w  -  /(e),  .. 
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3.3.2  Settlement  Success 


Behavioral  effects  on  settlement  success  can  be  considered  as  two  separate  questions:  1)  Un¬ 
der  a  given  set  of  physical  conditions,  which  behavior  is  most  successful  for  settlement?  and 
2)  For  a  given  behavior,  which  physical  conditions  allow  larvae  to  settle  most  successfully? 
Results  addressing  these  questions  are  discussed  separately  below. 

The  most  successful  behavior  for  given  conditions 

Larval  settlement  success  depends  on  combined  effects  of  behavior  and  the  timing  of  larval 
supply  and  attachment  probability  (Fig.  3-5).  Although  negative  buoyancy  always  results 
in  more  settlement  than  positive  buoyancy,  the  relative  benefits  of  turbulence-dependent 
behavior  depend  on  both  the  current  regime  and  the  settlement  function.  In  the  slowest 
currents  (Ur^  =  15  cm  s”1),  behaving  larvae  settle  less  successfully  than  passive  larvae 
(Fig.  3-5  d,  g)  because  they  are  more  positively  buoyant  and  have  lower  larval  supply 
throughout  the  tidal  cycle  (Figs.  3-1  e  and  3-4  a).  When  settlement  velocity  is  a  step 
function,  behaving  larvae  are  less  successful  than  positively-buoyant  larvae  at  all  current 
speeds  (Fig.  3-5  g-i)  because  behaving  larvae  are  swimming  upward  during  slack  tides  and 
have  very  low  larval  supply  when  settlement  occurs  (Figs.  3-3  and  3-4).  When  settlement 
velocity  is  linear  or  constant,  behaving  larvae  do  as  well  as  or  better  than  negatively-buoyant 
ones  at  moderate  (Ur^  =  35  cm  s'-1)  and  fast  {Ur^  =  65  cm  s^1)  current  speeds  (Fig. 
3-5  e,f).  The  relative  benefits  of  different  behaviors  are  consistent  with  depth  (but  see  Fig. 
3-5  e),  suggesting  that  these  results  can  be  generalized  for  larvae  settling  in  any  shallow 
habitats. 

The  most  favorable  conditions  for  a  given  behavior 

Passive  larvae  settle  most  successfully  in  the  slowest  currents  (Ur^  —  15  cm  s-1),  regard¬ 
less  of  whether  they  are  positively  or  negatively  buoyant  and  regardless  of  the  settlement 
function  (Fig.  3-6a-b).  It  is  intuitive  that  negatively-buoyant  larvae  would  settle  more 
successfully  when  the  diffusivity  is  low,  because  diffusivity  mixes  larvae  upward  that  would 
otherwise  sink  straight  to  the  bottom.  It  is  less  obvious  why  positively-buoyant  larvae  would 
settle  more  successfully  when  diffusivity  is  low,  because  diffusivity  mixes  larvae  downward 
that  would  otherwise  float  to  the  surface.  This  counter-intuitive  result  can  be  explained  by 
the  fact  that  larval  supply  is  similar  for  all  flow  speeds,  but  there  is  a  higher  time-averaged 
settlement  velocity  in  slower  flows. 

Unlike  passive  larvae,  larvae  with  turbulence-dependent  behavior  settle  most  successfully 
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Figure  3-6:  Proportion  of  larvae  settled  in  1  tidal  cycle  Q  vs.  current  velocity  U for 
three  settlement  velocity  functions  s.  Symbols  indicate  behavior  function:  w  =  0.05  cm 
s”1,  A;  w  =  —0.05  cm  s-1,  y;  w  =  /(e),  •. 

in  the  strongest  currents  (17//^  =  65  cm  s'"1)  for  the  constant  (not  shown)  and  linear 
settlement  functions  (Fig.  3-6c).  This  result  is  due  to  the  fact  that  larval  supply  is  high 
over  a  longer  time  period  as  current  speed  increases  (Fig.  3-4).  The  fact  that  behaving 
larvae  settle  more  successfully  in  faster  currents  is  interesting  because  the  time-averaged 
settlement  velocity  s  is  actually  lowest  in  the  strongest  currents.  This  result  indicates  that 
having  high  larval  supply  can  more  than  compensate  for  the  lower  settlement  velocities 
expected  in  more  turbulent  environments. 


44 


U  =  10  cm  s  1  U  =  65  cm  s 


Figure  3-7:  Larval  concentrations  from  PT  model  with  deterministic  behaviors  (y-axis)  vs. 
stochastic  behaviors  (x-axis)  for  two  flow  speeds.  Solutions  are  sampled  every  ten  minutes 
and  normalized  by  the  depth-integrated  concentration  Cj*.  Solid  lines  are  1:1. 
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Figure  3-8:  Larval  concentrations  from  PT  model  (y-axis)  vs.  AD  model  (x-axis)  for 
behaving  larvae  (w  =  /(e))  for  U =  65  cm  s_1  and  four  different  time-steps.  Solutions 
are  sampled  every  ten  minutes  and  normalized  by  the  depth-integrated  concentration  C^. 
Crosses  indicate  larval  concentrations  in  the  bottom  1  cm,  diamonds  indicate  concentrations 
in  the  upper  1  cm.  Solid  line  is  1:1. 
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Figure  3-9:  Larval  concentrations  from  PT  model  (y-axis)  vs.  AD  model  (x-axis)  for 
negatively  buoyant  larvae  with  (a)  w  =  —0.05  cm  s”1,  (b)  w  —  —0.2  cm  s~2,  and  (c) 
w  =  —0.5  cm  s_1  for  U =  65  cm  s'"1,  At  =  1.0  s.  Solutions  are  are  sampled  every  ten 
minutes  and  normalized  by  the  depth-integrated  concentration  Ct.  Crosses  indicate  larval 
concentrations  in  the  bottom  1  cm,  diamonds  indicate  concentrations  in  the  upper  1  cm. 
Dots  indicate  larval  concentrations  at  the  inner  grid  points  (1  <  z  <  100  cm)  Solid  lines  axe 
1:1. 
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3.3.3  Advection-Diffusion  vs.  Particle  Tracking  model 


The  particle  tracking  model  was  run  for  two  current  speeds  with  behaviors  specified  both 
implicitly  (3.10)  and  explicitly  (3.17),  and  the  larval  concentrations  for  the  two  model 
versions  are  very  similar  over  all  z  and  t  (Fig.  3-7).  Both  model  versions  contain  inaccuracies 
due  to  time-step  issues  (see  below),  but  the  solutions  are  affected  equally.  Given  the  nearly 
identical  results  for  implicit  and  explicit  behaviors,  it  is  safe  to  conclude  that  there  is  no 
loss  of  accuracy  when  behaviors  are  specified  implicitly  as  a  population-average  velocity. 

When  the  PT  solution  for  implicit  behavior  is  compared  against  the  AD  version  (Fig.  3- 
8),  it  is  clear  that  the  PT  solution  contains  significant  errors  near  the  boundaries  due  to  the 
time  step.  The  PT  solution  deviates  less  from  the  AD  solution  as  the  time-step  is  reduced, 
and  presumably  would  be  accurate  if  the  time-step  were  small  enough.  The  PT  model 
(3.15)  is  equivalent  to  the  AD  model  (3.1)  only  as  the  number  of  particles  N  — >  oo  and  as 
At  — ♦  0  [7].  For  At  0,  the  stochastic  terms  in  (3.1o)  ensure  that  particles  will  sometimes 
travel  a  long  distance  and  bounce  off  the  boundaries.  The  result  is  an  underestimate  of 
concentration  at  the  boundaries,  as  seen  in  Figure  3-8. 

The  underestimation  of  near-boundary  concentration  is  more  pronounced  as  the  advec- 
tive  velocity  w  grows  larger  in  magnitude  (Fig.  3-9).  When  both  models  were  run  for 
passive  larvae  with  three  different  values  of  w  (=  [-0.05,  -0.2,  -0.5]  cm  s-1),  the  PT  solution 
was  reasonably  accurate  only  when  larvae  were  nearly  neutrally-buoyant  (w  —  —0.05). 

For  the  hydrodynamic  conditions  and  behaviors  discussed  in  this  paper,  the  AD  model 
is  preferable  to  the  PT  model  for  two  reasons.  First,  behavioral  velocities  are  small  enough 
relative  to  the  diffusivity  that  multiple  behaviors  can  be  modeled  implicitly  as  a  population- 
average  velocity  with  no  loss  of  accuracy.  Second,  the  PT  model  requires  very  small  time 
steps  to  minimize  errors,  particularly  for  larger  w .  Although  the  PT  model  is  attractive  in 
allowing  greater  flexibility  in  modeling  larval  behavior,  it  cannot  be  counted  on  to  produce 
accurate  results  near  the  boundaries  over  the  short  spatial  scales  modeled  here. 


3.4  Discussion 
3.4.1  Larval  Supply 

Turbulence-induced  sinking  behavior  significantly  affects  the  temporal  patterns  of  larval 
supply  to  the  bed.  In  moderate  to  strong  currents,  negatively-buoyant  larvae  have  large 
peaks  in  near-bed  concentration  during  slack  tides,  whereas  behaving  larvae  are  highly 
concentrated  at  the  bed  during  flood  and  ebb  tides.  These  opposite  patterns  of  larval 
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supply  could  each  be  optimal  for  settlement  under  different  environmental  conditions. 

Intuitively  it  might  seem  that  larvae  would  benefit  most  from  concentrating  near  the 
bottom  during  slack  tides,  when  shear  stress  is  low  and  sediments  are  stable.  Yet  as  pointed 
out  by  Eckman  et  al  [8],  there  is  a  cost  to  reaching  the  bottom  during  slack  tide  if  the 
substrate  is  unsuitable  for  settlement.  Larvae  that  reject  substrates  at  slack  tide  have 
to  wait  for  currents  to  increase  again  and  carry  them  away  to  potentially  better  sites. 
Furthermore,  larvae  that  contact  the  bottom  only  during  slack  tides  might  have  just  one 
settlement  opportunity  per  tidal  cycle.  In  contrast,  larvae  that  reach  the  bottom  during 
flood  and  ebb  tides  could  benefit  from  being  able  to  test  substrates  and  still  be  transported 
rapidly  away  from  unsuitable  sites.  Larvae  with  a  high  near-bed  concentration  during  most 
of  the  tidal  cycle  could  also  have  more  settlement  opportunities,  for  example  during  lulls  in 
turbulence  (discussed  below).  Negative  buoyancy  results  in  peaks  of  larval  supply  at  slack 
tides,  and  might  confer  the  optimal  timing  of  larval  supply  in  a  world  where  all  substrates  are 
suitable  for  settlement.  Turbulence-induced  sinking  results  in  high  larval  supply  throughout 
flood/ebb  tides,  and  is  potentially  a  more  successful  strategy  in  patchy  environments. 

3.4.2  Settlement  Success 

In  the  simplified,  1-dimensional  environment  explored  in  this  paper,  larvae  that  sink  in  tur¬ 
bulence  are  more  successful  settlers  than  passive  larvae  under  some  conditions.  In  moderate 
to  fast  currents  {JJhoo  ~  35  and  65  cm  s-1)  turbulence-induced  sinking  is  the  most  success¬ 
ful  behavior  for  settlement  as  long  as  the  settlement  velocity  s  is  non-zero  during  flood/ebb 
tides  (constant  and  linear  s).  When  the  settlement  velocity  is  zero  during  peaks  of  larval 
supply  (during  flood  and  ebb  tides,  step  s),  then  almost  no  settlement  occurs.  The  tim¬ 
ing  of  near-bed  hydrodynamics  is  therefore  key  to  determining  whether  turbulence-induced 
sinking  can  be  a  more  successful  settlement  strategy  than  passive  transport. 

The  simplified  settlement  functions  in  this  model  ignored  two  important  boundary  layer 
characteristics  that  would  probably  improve  the  relative  settlement  success  of  behaving 
larvae.  First,  the  model  excluded  turbulence  intermittency  at  the  bed.  Even  during  flood 
and  ebb  tides,  there  are  lulls  between  turbulent  bursts  at  the  bed.  The  duration  of  these 
lulls  can  be  estimated  as  ~  6 H/Uh^  [5,  22],  and  at  peak  flood/ebb  tides  would  be  on  the 
order  of  ~  10  to  70  s  for  the  depths  and  current  speeds  used  in  this  study.  The  shear 
stress  is  lower  during  lulls,  and  larvae  would  almost  certainly  be  able  to  stick  or  burrow  into 
sediments  during  some  of  these  periods.  Behaving  larvae  that  have  high  larval  supply  during 
flood/ebb  tides  are  probably  able  to  exploit  the  intermittent  stress  lulls  for  settlement. 
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The  second  boundary-layer  simplification  in  this  model  was  the  absence  of  roughness 
elements.  Roughness  elements  influence  the  vertical  diffusivity  profile  near  the  bed  as  well 
as  the  horizontal  distribution  of  shear  stress  along  the  bottom.  Although  the  roughness- 
element  spacing  has  complex  effects  on  larval  attachment  probability  [5],  larval  settlement 
is  generally  greater  over  dense  roughness  elements  than  over  flat  beds  [9].  Flat  beds  are 
rare  in  tidal  inlets  such  as  Barnstable  Harbor,  where  ripples  and  epifauna  provide  small- 
scale  bottom  topography.  Mud  snails  themselves  form  dense  aggregations  with  a  roughness 
height  of  1-2  cm  (Fig.  1-1),  and  their  presence  could  alter  the  diffusivity  and  shear  stress 
profiles  to  enhance  settlement  of  larvae  during  flood/ebb  tides. 

Modeled  larvae  were  allowed  no  behavioral  responses  to  environmental  variables  other 
than  turbulence,  but  this  is  clearly  an  oversimplification.  Mud  snail  larvae  settle  selectively 
in  sediment  from  adult  habitats  [27],  and  changes  in  behavior  near  the  bed  could  increase 
both  the  supply  of  larvae  and  the  attachment  probability.  Larvae  within  a  short  distance 
(e.g.  <  1  cm)  of  the  bottom  probably  sense  biological  and  chemical  properties  of  the  sedi¬ 
ments,  and  might  have  a  lower  probability  of  swimming  upward  when  these  sediments  are 
attractive  for  settlement.  This  would  result  in  a  more  negative  average  larval  velocity  near 
the  bed,  and  larvae  would  tend  to  pile  up  near  the  bottom.  In  the  laboratory,  larvae  settle 
by  burrowing  into  sediments;  this  burrowing  behavior  could  increase  the  settlement  veloc¬ 
ity  by  increasing  the  shear  velocity  at  which  larvae  are  eroded.  Any  additional  behavioral 
responses  would  be  likely  to  increase  larval  settlement  success,  although  it  remains  unclear 
whether  larvae  that  sink  in  turbulence  would  have  greater  overall  settlement  success  than 
passive  larvae  if  substrate  selection  were  allowed. 

Finally,  this  model  does  allow  some  general  predictions  about  where  larvae  are  most 
likely  to  settle,  given  their  behavior.  Although  negatively  buoyant  larvae  settle  successfully 
under  all  conditions,  the  settlement  of  passive  larvae  is  most  likely  to  occur  in  slow  to 
moderate  currents.  Thus  for  intertidal  species  that  prefer  calm  or  sheltered  habitats,  no 
response  to  turbulence  is  necessary.  On  the  other  hand,  larvae  with  turbulence-dependent 
behavior  are  more  likely  to  settle  in  faster  currents  such  as  those  in  Barnstable  Harbor. 
Mud  snails  are  abundant  in  Barnstable  Harbor,  and  it  is  tempting  to  conclude  that  mud 
snail  larvae  sink  in  turbulence  in  order  to  settle  into  preferred  habitats.  However  it  is  still 
impossible  to  reject  the  alternative  hypothesis  that  gastropod  larvae  sink  in  turbulence  to 
avoid  predators..  Regardless  of  the  ultimate  reason  for  sinking,  it  is  clear  that  larvae  with 
this  behavior  could  experience  enhanced  settlement  in  turbulent  coastal  inlets  relative  to 
quieter  inlets  or  calmer  offshore  waters. 
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Chapter  4 


Larval  behavior  in  turbulence: 
Maximum-likelihood  estimates 
from  field  plankton  distributions 


4.1  Introduction 

Larval  behavior  contributes  to  larval  vertical  distributions,  which  in  turn  influence  dispersal 
distances  and  the  supply  of  larvae  to  benthic  habitats.  Larval  supply  is  correlated  to 
settlement  at  least  over  ~  weekly  time  scales  [24,  37],  and  behavioral  effects  on  larval  supply 
are  important  insofar  as  larval  supply  can  be  considered  a  proxy  for  settlement.  Behaviors 
that  influence  the  delivery  of  larvae  to  benthic  habitats  are  still  poorly  understood,  but  are 
potentially  very  important  in  determining  spatial  and  temporal  settlement  patterns. 

Most  of  our  knowledge  about  larval  behavior  comes  from  laboratory  studies,  but  these 
studies  provide  limited  insight  into  larval  behavior  under  natural  conditions.  Turbulence  [9, 
27,  45]  and  many  other  physical  [e.g.,  8,  10,  39]  and  chemical  [e.g.,  17,  26,  40]  cues  elicit  larval 
behaviors  in  the  laboratory.  In  the  natural  environment  however,  larvae  encounter  multiple, 
simultaneous  cues  that  change  over  different  spatial  and  temporal  scales.  Taken  together, 
these  multiple  cues  could  suppress  or  enhance  behaviors.  Turbulence  has  been  identified  as 
a  potentially  important,  settlement-related  cue  for  blue  crabs  [45]  and  mud  snails  ( Ilyanassa 
obsoleta )  [9]  in  the  laboratory,  but  the  effect  of  turbulence  on  larval  behavior  has  yet  to  be 
tested  in  the  environment. 
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Estuaries  are  an  ideal  place  to  field-test  whether  larvae  respond  to  turbulence.  Tur¬ 
bulence  in  estuaries  changes  rapidly  over  tidal  periods  and  varies  over  short  spatial  scales 
(in).  Most  other  physical  cues  (e.g.,  pressure,  temperature,  salinity)  also  vary  over  rela¬ 
tively  short  temporal  (h)  and  spatial  scales  (m  to  km),  but  none  vary  in  space  and  time 
in  the  same  way  turbulence  does.  Therefore  responses  to  other  physical  cues  are  unlikely 
to  be  mistaken  for  a  response  to  turbulence.  Larval  behaviors  may  depend  on  multiple 
physical  factors,  all  of  which  vary  in  estuaries.  If  a  larval  response  to  turbulence  is  detected 
in  this  dynamic  environment,  then  turbulence  could  be  taken  as  an  important  determinant 
of  overall  behavior. 

The  present  study  was  conducted  at  Barnstable  Harbor,  MA  (USA),  in  order  to  estimate 
larval  responses  to  turbulence  in  the  field.  Barnstable  Harbor  is  a  well- mixed  coastal  inlet 
with  low  freshwater  inflow.  The  boundary  layer  is  fully  turbulent  and  depth-limited,  with 
maximum  current  speeds  of  50  to  80  cm  s-1  [1,  20].  The  tidal  range  averages  ~3  m,  and  60% 
of  the  bottom  is  exposed  at  low  tide.  The  intertidal  flats  are  covered  with  well-sorted,  fine 
sand  and  silt,  with  an  average  sediment  grain  diameter  of  ~0.1  mm  throughout  the  harbor 
[35,  33].  Mud  snails  are  the  most  conspicuously  abundant  megafauna  on  the  intertidal  flats, 
and  were  a  focus  of  this  study. 

We  formed  two  hypotheses  about  mud  snail  larval  behavior  in  the  field,  based  on  our 
previous  work  on  mud  snail  larvae.  In  the  laboratory,  competent  (able  to  metamorphose) 
mud  snail  larvae  exhibit  three  behaviors:  swimming,  hovering,  and  sinking  with  the  velums 
withdrawn.  The  proportion  of  sinking  larvae  increases  exponentially  with  the  turbulence 
dissipation  rate  £,  and  on  average  larvae  switch  from  upward-swimming  to  sinking  at  dissi¬ 
pation  rate  of  about  10_1  cm2  s~3.  Modeling  results  suggest  that  this  behavior  will  enhance 
settlement  in  coastal  inlets  like  Barnstable  Harbor.  We  hypothesize  that  competent  larvae 
have  a  similar  response  to  turbulence  in  the  field,  and  sink  more  in  stronger  turbulence.  Be¬ 
cause  sinking  in  turbulence  is  thought  to  be  a  settlement-related  behavior,  we  hypothesize 
that  pre-competent  larvae  have  little  or  no  response  to  turbulence. 

We  also  used  this  study  as  an  opportunity  to  ask  the  question,  do  larvae  respond  dif¬ 
ferently  to  turbulence  on  flood  and  ebb  tides?  Tidally-dependent  behaviors  allow  larvae  of 
some  estuarine  crab  species  to  exit  or  enter  coastal  embayments  at  appropriate  life  stages 
[5,  45].  Mud  snail  larvae  might  have  different  settlement  success  if  they  sink  to  the  bottom 
during  flood  vs.  ebb  tides.  Although  our  laboratory  results  allow  us  to  predict  that  mud 
snail  larvae  will  sink  more  in  stronger  turbulence,  our  experiments  were  conducted  at  con¬ 
stant  temperature  and  salinity,  giving  us  no  empirical  basis  to  predict  whether  this  behavior 
depends  on  tidal  stage.  Our  field  study  was  done  in  a  tidal  environment,  allowing  us  to  test 
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whether  tidal  stage  affects  the  larval  response  to  turbulence. 

Larval  behaviors  are  difficult  to  quantify  in  controlled  experiments,  let  alone  in  uncon¬ 
trolled  natural  habitats.  Larval  ecologists  often  use  the  vertical  distributions  of  larvae  in 
the  field  to  speculate  about  larval  behavior  [e.g.,  3,  5,  29].  Remote  sampling  techniques 
provide  more  quantitative  ways  to  observe  and  estimate  plankton  velocities  in  the  field 
[11,  12],  but  species-level  analyses  of  remotely-sampled  data  are  still  difficult  or  impossi¬ 
ble.  In  this  study,  behaviors  were  estimated  quantitatively  from  the  vertical  distributions 
of  larvae  in  a  turbulent  boundary  layer,  rather  than  from  behavioral  observations.  Larval 
velocities  were  estimated  as  a  function  of  the  turbulence  dissipation  rate,  using  a  maximum 
likelihood  analysis  of  larval  concentration  distributions.  These  maximum-likelihood  behav¬ 
ior  estimates  represent  population-level  responses  to  turbulence.  The  results  suggest  that 
other  cues  are  also  important,  but  turbulence  does  elicit  a  strong  behavioral  response  for 
crab  and  gastropod  larvae  in  the  field. 


4.2  Methods 

4.2.1  Plankton  Samples 

Plankton  samples  and  CTD  data  were  collected  from  a  small  boat  moored  in  the  Barnstable 
Harbor  channel  (Fig.  4-1)  in  summer  2004  (22,  27,  29  July).  Each  day  involved  10  horns  of 
sampling  beginning  at  07:30.  Hourly  measurements  consisted  of  one  CTD  profile  and  one 
vertical  profile  of  plankton  samples.  CTD  profiles  were  collected  with  an  Ocean  Sensors 
OS200.  Plankton  samples  were  collected  at  the  surface  and  in  1-m  increments  to  the 
bottom  (Fig.  4-2),  using  an  Ebara  dominator  pump  with  a  Great  Plains  Industries  TM200 
flowmeter.  A  Scotty  downrigger  was  used  to  set  the  plankton  pump  at  precise  depths.  For 
each  plankton  sample,  300  L  of  water  was  pumped  through  a  200  /zm-mesh  plankton  net. 

Plankton  samples  were  stored  overnight  at  4  °C  and  preserved  on  the  following  day 
except  the  22  July  samples,  which  were  preserved  2  days  later.  Samples  were  fixed  in  4% 
formalin  with  borax  buffer  for  up  to  2  weeks,  and  then  transferred  to  buffered  80%  ethanol. 
For  sorting,  22  July  samples  were  sieved  into  200-500  n m  and  >500  fim  size  fractions. 
All  other  samples  were  sieved  into  200-425  ^m  and  >425  jj. m  size  fractions.  These  size 
fractions  correspond  to  the  approximate  size  ranges  of  pre-competent  and  competent  mud 
snail  larvae.  Samples  from  22  July  were  poorly  preserved  and  were  not  identified  to  species. 

Invertebrate  larvae  were  sorted  and  identified  to  major  taxa.  Only  brachyuran  crab 
zoea  and  gastropod  larvae  were  found  in  most  samples  and  were  counted.  Because  we 
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Figure  4-1:  Cape  Cod  map  showing  location  of  Barnstable  Harbor,  and  Barnstable  Harbor 
bathymetry  map  with  locations  of  plankton  mooring  (PM)  and  current  meters  (CM).  Dark 
outlines  indicate  mean  high  water  and  white  lines  indicate  mean  low  water  (MLW). 

were  interested  in  behaviors  specific  to  gastropod  larvae,  gastropods  were  identified  to 
species  using  descriptions  in  the  literature  [34,  36,  41].  In  addition  to  Ilyanassa  obsoleta 
larvae,  Crepidula  spp.  (including  Crepidula  plana  and  Crepidula  fomicata)  and  Anachis 
spp.  (including  Anachis  avara  and  Anachis  translirata)  were  found  in  at  least  half  the 
profiles  and  were  used  for  a  genus-level  analysis.  Other  species  of  gastropods  (e.g.,  Bittium 
altematum ,  Lunatia  heros)  were  too  rare  for  a  species-  or  genus-level  behavioral  analysis. 

4.2.2  Current  Measurements 

Flow  velocities  were  measured  with  an  acoustic  Doppler  current  profiler  (ADCP)  and  an 
acoustic  Doppler  velocimeter  (ADV),  deployed  in  the  Barnstable  Harbor  channel  (Fig.  4-1) 
during  the  plankton  sampling  period.  The  ADCP  measured  current  profiles  but  was  unable 
to  measure  velocities  within  about  1.5  m  of  the  bottom,  and  the  ADV  provided  an  additional 
measurement  closer  to  the  bottom.  The  instruments  were  deployed  20  July  and  recovered 
17  August  2004  from  the  R/V  Tioga.  A  1.2-MHz,  RDInstruments  Workhorse  ADCP  was 
mounted  on  a  tripod,  with  the  transducer  heads  looking  up  and  located  ~0.75  m  above  the 
bed.  The  tripod  was  leveled  to  a  ~2°  tilt,  wuthin  the  acceptable  limits  for  Reynolds  stress 
calculations  [22].  Approximately  10  m  south  of  the  ADCP,  a  5-MHz,  Sontek  ADVOcean 
was  mounted  by  divers  to  a  pipe  jetted  into  the  sand,  with  the  transducer  head  pointed 
up  and  located  ~0.5  m  above  the  bed.  The  channel  bottom  was  covered  with  large  sand 


58 


waves  (height  ~  2  m,  wavelength  ~  6  to  10  m)  aligned  perpendicular  to  the  channel.  Both 
instruments  were  deployed  atop  the  ridge  of  a  sand  wave. 

The  ADCP  collected  data  continuously  from  about  15:00  20  July  to  about  12:00  14  Au¬ 
gust.  Using  RDI’s  mode  12,  10-Hz  sub-pings  were  averaged  over  one  second  and  recorded 
every  other  second.  Velocities  were  recorded  in  beam  coordinates  for  Reynolds  stress  esti¬ 
mation.  To  resolve  the  shallow  water  column,  a  25-cm  bin  size  was  used,  with  the  first  bin 
centered  at  1.5  m  from  the  bed.  The  ADV  collected  data  from  00:00  22  July  until  14:00  17 
August.  Velocities  were  measured  in  earth  coordinates  ( E ,  N ,  and  Up)  at  8  Hz  for  512  s 
bursts  every  half  hour.  The  ADV  measurement  volume  was  located  about  0.7  m  from  the 
bed.  ADCP  and  ADV  data  are  reported  for  22  July  to  31  July. 

4.3  Analysis 

4.3.1  Current  Measurements 
ADCP 

The  ADCP  was  used  to  measure  current  and  turbulence  profiles  in  the  water  column.  For 
the  calculation  of  mean  current  velocities  and  shear,  velocities  were  transformed  from  beam 
to  cartesian  coordinates  and  corrected  for  pitch  and  roll  [21,  38].  Currents  were  rotated  to 
coordinates  where  u(z,t)  is  the  along-channel  velocity  (positive  towards  Cape  Cod  Bay), 
v(z,t)  is  the  cross-channel  velocity,  and  w(z,t)  is  the  vertical  velocity  (positive  upwards). 
Velocities  for  each  depth  bin  were  calculated  as  10-minute  averages,  U,  V ,  and  W.  The 
depth-averaged  shear  velocity  u+d  and  hydraulic  roughness  length  z0  were  calculated  from 
the  along-current  velocity  profiles  using  the  Law  of  the  Wall  [e.g.,  4], 

=  (4.1) 

&  '\^o  / 

where  k  is  von  Karman’s  constant  (=  0.4).  A  depth- averaged  drag  coefficient  was  calculated 
as  Cd  =  ulJU2,  where  U  is  the  depth- averaged  along-current  velocity  [16]. 

Reynolds  stress  u'w'  was  calculated  from  the  original  ADCP  beam  velocities  over  10- 
minute  intervals  using  the  variance  method  [e.g.,  22,  38].  Relevant  turbulence  quantities 
were  calculated  as  functions  of  the  Reynolds  stress  and  the  along-current  velocity  shear. 
The  vertical  eddy  diffusivity  Kp  is  defined  as 

Kp  =  ~w/^.  (4.2) 
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The  turbulence  dissipation  rate  ev  was  assumed  to  be  equal  to  the  production  of  turbulent 
kinetic  energy  [e.g.,  43],  estimated  as 


£  =  — utu 


7dU 

Hz' 


(4.3) 


ADV 

The  ADV  was  used  to  measure  3-dimensional  velocities  near  the  bed.  The  velocity  data 
included  bad  measurements  due  to  obstruction  of  the  sensor  head  by  crabs  or  floating  algae 
[P.  Schultz,  pers.  comm.].  The  data  were  pre-processed  to  remove  low-quality  measurements 
[e.g.,  7],  and  29%  of  recorded  bursts  were  thrown  out.  Velocities  were  then  rotated  to 
coordinates  where  u(t)  is  the  along-channel  velocity,  v(t)  is  the  cross-channel  velocity,  and 
w(t)  is  the  vertical  velocity. 

The  mean  velocities,  shear  velocity,  drag  coefficient,  and  turbulence  dissipation  rate 
were  calculated  for  each  512-s  burst.  Velocity  time  series  were  de-trended  and  de-meaned 
for  calculation  of  turbulence  quantities.  The  shear  velocity  n*  was  calculated  directly  from 
the  Reynolds  stress  as  it*  =  |u/u;/|0‘5,  where  u'  and  wf  are  the  along-current  and  vertical 
velocity  fluctuations.  A  drag  coefficient  was  calculated  for  each  burst  as  =  u2/f/2, 
where  U  is  the  average  velocity  [16].  Dissipation  rate  estimates  ew  were  obtained  from  the 
one-sided  wavenumber  spectra  Ew  of  vertical  velocity  by  fitting 

B.  =  a  (|)  cT  (*-5/3  +  <*.  -  *)-s/3)  (4.4) 

in  the  inertial  subrange  [43],  where  a  ~  1.5  is  an  empirical  constant  [14],  k  is  the  wavenumber 
(related  to  frequency  /  by  k  =  27t//|(/|),  and  ks  is  the  wavenumber  corresponding  to  the 
sampling  frequency.  The  w  spectra  were  unaffected  by  noise  but  showed  some  effects  of 
aliasing  at  higher  flow  velocities,  and  the  wavenumber  term  of  (4.4)  includes  a  correction 
for  aliasing. 


4.3.2  Advection-Diffusion  Model 

Larval  behaviors  were  quantified  in  terms  of  vertical  velocities;  although  behaviors  can  be 
characterized  in  other  ways,  it  is  the  swimming  or  sinking  velocities  that  determine  vertical 
advection  of  larvae  toward  or  away  from  settlement  sites.  Larval  vertical  velocities  were 
estimated  as  a  function  of  the  turbulence  dissipation  rate  by  fitting  a  vertical  advection- 
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diffusion  (AD)  model 


dc  _  d  (  ec 

dt  dz  Ks  dz 


(4.5) 


to  the  observed  plankton  distributions.  This  model  predicts  the  concentration  of  larvae 
C(Zjt)  at  height  z  and  time  t.  Larvae  are  advected  due  to  their  behavioral  velocity  w(z,t) 
and  diffused  due  to  turbulent  eddy  diffusivity  K3(z,t).  Boundary  conditions  are  no-flux, 
so  larvae  bounce  off  the  bottom  and  surface.  No  settlement  was  included  in  the  boundary 
conditions  because  the  mixing  time  scale  was  assumed  to  be  shorter  than  the  settlement 
time  scale. 

The  larval  velocity  was  modeled  as  a  linear  function  of  the  turbulence  dissipation  rate, 


w  =  60  +  bx  \og10£s. 


(4.6) 


In  this  model  bi  is  the  most  interesting  parameter,  as  it  predicts  whether  larval  velocity 
is  positively  or  negatively  related  to  the  dissipation  rate.  If  larvae  have  no  response  to 
turbulence,  then  b\  —  Os2  cm-1  and  bo  indicates  whether  larvae  are  positively  or  negatively 
buoyant. 

In  order  to  fit  the  AD  model  to  the  plankton  data,  (4.5)  had  to  be  run  with  diffusivities 
Ks  and  dissipation  rates  £s  that  were  accurate  for  each  plankton  sampling  period.  The 
numerical  solution  of  (4.5)  also  required  that  Ks  and  £s  be  continuous  in  time.  There¬ 
fore,  high-order  polynomial  functions  were  fitted  to  field  data  to  produce  diffusivity  and 
dissipation  functions  that  were  both  accurate  and  continuous  (Appx.  B). 

The  AD  model  was  run  for  a  range  of  parameters  (—0.3  <  bo  <  0.2  cm  s-"1;  —0.3  <  b\  < 
0.3  s2  cm^1)  that  were  reasonable  given  the  range  of  larval  velocities  in  the  literature  [e.g. 
2,  9]  and  the  expected  range  of  dissipation  rates  [e.g.,  13,  16,  43].  A  water  depth  of  H  =  4 
m  was  used  in  the  model,  approximating  the  mean  depth  at  the  plankton  sampling  site. 
A  constant-depth  model  closely  approximated  the  depth-normalized  distribution  of  larvae 
in  a  variable-depth  water  column  (see  Sec.  4.3.4).  The  model  was  run  separately  for  each 
plankton  sampling  day,  using  the  appropriate  diffusivity  and  dissipation  functions.  Results 
were  used  in  the  maximum  likelihood  analysis  described  in  the  next  section. 


4.3.3  Maximum  Likelihood  Estimation  of  Behavior 

Plankton  data  were  analyzed  separately  for  each  size  class  and  5  groups  of  larvae:  Ilyanassa 
obsolete  Crepidula  spp.,  Anachis  spp.,  pooled  gastropod  larvae,  and  pooled  Brachyuran 
crab  zoea.  Plankton  profiles  from  flood  and  ebb  tides  were  analyzed  separately  (Fig.  4-2), 
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to  test  for  a  dependence  of  behavior  on  tidal  stage. 


The  AD  model  (4.5)  was  fitted  to  the  plankton  data  by  maximum  likelihood.  The 
likelihood  Lj{6)  of  each  observed  plankton  profile  was  calculated  for  each  combination  of 
parameter  values  6  =  ( bo,b\ )  using  the  multinomial  model 

=  Njl  biiWnjl^2(g)n^...Pj/c(g)n^] 

where  Nj  is  the  total  number  of  larvae  in  profile  j,  rij ^  is  the  number  of  larvae  in  sample 
k  of  profile  j  of  the  data,  and  Pjk{9)  is  the  proportion  of  larvae  in  sample  k  of  profile  j  of 
the  model,  given  6.  The  Pjk{&)  were  found  by  samphng  the  AD  model  solution  at  the  same 
proportional  depths  Zjz/H)  time,  and  date  corresponding  to  plankton  profile  j,  and  dividing 
by  the  total  concentration  in  the  model  profile.  Larval  behavior  parameters  6  —  (&o,&i) 
were  estimated  as  those  that  maximized  the  log  likelihood  LL(6) 

J 

LL(9)  -  ElniiW  (4-8) 

l 

over  all  profiles  in  each  group.  A  95%  confidence  interval  for  the  slope  was  estimated 
using  the  profile  likelihood. 


A  likelihood  ratio  was  used  to  test  the  hypothesis  that  the  larval  response  to  turbulence 
is  different  for  flood  and  ebb  tides, 


Ho  ■  $flood  —  $ebb  ,  . 

(4.9) 

H A  :  $flood  7*”  $ebb 

where  0flood  and  $ebb  are  the  behavior  parameters  for  flood  and  ebb  tides.  The  null  hypoth- 
esis  was  rejected  if  2 [LL(6Rood)  +  LL(6ehb)  -  LL{dboo d+ebb)]  >  xl.0.05)  where  0flood  and  6ebb 
are  maximum-likelihood  parameter  estimates  (MLE’s)  for  flood-tide  and  ebb-tide  samples, 
and  0fiood+ebb  are  the  MLEs  for  combined  flood-  and  ebb-tide  samples.  This  hypothesis 
was  tested  separately  for  the  small  and  large  size  fractions  of  each  larval  group.  A  similar 
test  was  used  for  the  hypothesis  that  large  and  small  size-fraction  larvae  respond  differently 
to  turbulence.  A  Chi-square  goodness-of-fit  test  was  used  to  determine  whether  the  linear 
behavior  model  (4.6)  was  correct  for  each  larval  group,  tidal  stage,  and  size  fraction. 
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Figure  4-2:  Schematic  of  plankton  samples  collected  on  (a)  22  July,  (b)  27  July,  and  (c) 
29  July.  Symbols  show  classification  of  samples  as  flood  or  ebb  tide  profiles  or  transitional 
profiles  (analysed  with  both  flood  and  ebb  tide  samples).  Two  profiles  were  discarded  from 
the  analysis  due  to  stratification  (>1  °C  difference  between  surface  and  bottom).  Thick 
solid  line  indicates  water  depth  relative  to  the  surface. 
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4.3.4  Model  assumptions 

Behavior  estimates  were  made  from  the  vertical  distributions  of  larvae  in  a  turbulent  bound¬ 
ary  layer,  rather  than  from  behavioral  observations.  This  analysis  required  several  assump¬ 
tions  about  how  larvae  are  transported  in  turbulent  eddies. 

A  key  assumption  of  this  analysis  is  that  neutrally  buoyant  particles  are  evenly  dis¬ 
tributed,  and  any  deviations  from  an  even  vertical  distribution  are  due  to  behavior.  This 
assumption  is  valid  only  if  the  water  column  is  well  mixed,  as  stratification  could  indicate 
the  presence  of  separate  water  masses  containing  isolated  larval  patches  [e.g.,  10].  Plankton 
profiles  were  excluded  from  the  analysis  if  the  concurrent  CTD  profiles  showed  any  strat¬ 
ification,  defined  as  a  temperature  difference  of  >  1  °C  between  the  surface  and  bottom. 
Barnstable  Harbor  is  very  well-mixed,  and  by  this  criteria  only  two  plankton  profiles  had 
to  be  excluded. 

Another  assumption  is  that  larval  diffusivity  is  equivalent  to  fluid  diffusivity.  This 
assumption  is  valid  if  the  transport  of  larvae  by  eddies  is  exempt  from  particle  inertia  and 
crossing-trajectory  effects  [e.g.,  44,  32].  These  effects  were  considered  negligible  based  on 
criteria  reviewed  in  Ross  and  Sharpies  [32].  Large,  dense  particles  have  inertia  that  increases 
the  response  time  rp  or  lag  time  between  eddy  motion  and  motion  of  a  particle  [32,  Eq.  10]. 
The  response  time  of  a  1  mm  larva  with  a  density  of  1.07  g  cm”3  is  rp  <  0.1  s,  and  inertial 
effects  are  negligible.  Trajectory-crossing  occurs  when  particles  have  vertical  velocities  wp 
greater  than  the  eddy  velocity  fluctuations  w' .  Fast  particles  fall  through  eddies  rather 
than  following  them,  and  particle  motion  becomes  uncorrelated  with  the  fluid  motion.  The 
correlation  of  particle  motion  and  fluid  motion  can  be  quantified  as  the  ratio  of  particle 
diffusivity  Kp  to  fluid  diffusivity  Kj  [32,  Eq.  11].  For  an  average  larva  (wp  =  0.5  cm  s”1) 
in  Barnstable  Harbor,  Kp/Kj  >  0.97  and  crossing-trajectory  effects  are  negligible. 

It  is  also  assumed  that  the  AD  model  with  constant  depth  H  accurately  predicts  the 
particle  distribution  in  a  water  column  of  variable  depth.  This  assumption  was  tested  using 
a  particle-tracking  equivalent  of  the  advection  diffusion  model  (Chapter  3).  The  model  was 
run  for  three  behavior  functions  (w  =  —0.05  cm  s”1,  w  =  0.05  cm  s”1,  w  =  /(e))  with 
constant  depth  (Hc  =  4  m)  and  with  a  1.5  m  tidal  amplitude  (Hv(t)  —  4  m)  over  a  full  tidal 
cycle.  Results  were  sampled  at  the  same  proportional  depths.  A  least-squares  regression 
confirmed  that  a  constant-depth  model  accurately  predicts  the  distribution  of  particles  in 
water  of  variable  depth  ( R 2  =  0.99). 

Finally,  this  analysis  required  the  assumption  that  Barnstable  Harbor  currents  are  later¬ 
ally  homogeneous  over  small  (100’s  m)  spatial  scales,  because  the  plankton  data  and  hydro 
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dynamic  data  were  collected  about  500  m  apart.  An  assumption  of  laterally-homogeneous 
flow  is  used  frequently  in  plankton  studies  [e.g.,  12]  because  a  separation  between  instru¬ 
ments  is  necessary  to  avoid  flow  disturbance.  In  this  study  the  sites  were  selected  in  co¬ 
operation  with  the  harbor  master  to  avoid  navigational  hazards,  and  the  distance  between 
sites  was  unavoidable.  A  potential  problem  with  the  lateral  homogeneity  assumption  is 
addressed  in  the  discussion. 


4.4  Results 

4.4.1  Current  Measurements 

The  water  column  at  Barnstable  Harbor  was  well  mixed,  as  evidenced  by  nearly  vertical 
temperature  and  salinity  profiles  (Fig.  4-3).  Salinity  varied  by  no  more  than  about  1 
psu  over  any  10-hour  sampling  period,  indicating  low  freshwater  inflow  into  the  harbor. 
Temperature  varied  by  5  to  6  °C  between  high  tide  and  low  tide,  due  to  rapid  warming  of 
the  shallow  inlet  during  low  tides. 

Maximum  current  velocities  were  about  80  cm  s'-1  when  averaged  over  10  min  intervals 
(Fig.  4- 4 a),  with  higher  velocities  on  flood  tides  than  on  ebb  tides.  Shear  velocity  estimates 
from  the  ADCP  velocity  profiles  were  up  to  10  cm  s_1  (Fig.  4-4b),  whereas  the  maximum 
estimate  from  the  ADV  was  only  6.6  cm  s-1  (Fig.  4-5b).  It  is  likely  that  the  velocity  profile 
method  overestimated  the  shear  velocity.  Drag  coefficients  from  the  ADV  measurements 
averaged  C d  =  1.34  x  10~2,  with  higher  values  on  ebb  tides  than  on  flood  tides  (Fig.  4-5b). 

Dissipation  rates  were  tidally  asymmetric,  and  high  relative  to  values  reported  for  other 
tidal  inlets  (Fig.  4-5d).  Near  the  bed,  dissipation  rates  €w  from  ADV  velocity  spectra 
ranged  from  10-3  2  to  101  0  cm2  s-3  with  higher  values  on  ebb  tides  than  on  flood  tides 
(Fig.  4-5d).  The  tidal  signal  was  less  apparent  in  surface  estimates  of  dissipation  rate  ep 
(Fig.  4-4d),  which  were  about  2  orders  of  magnitude  lower  than  near-bottom  estimates. 
Tidal  asymmetries  were  captured  well  by  the  smooth  function  es  that  was  fitted  to  ep  and 
used  in  the  larval  behavior  analysis  (Fig.  4-7). 

Eddy  diffusivity  estimates  Kp  were  also  tidally  asymmetric  and  generally  high  (Fig.  4- 
6).  Flood  tide  diffusivity  profiles  had  typical  mid-depth  maxima  of  a  few  hundred  cm2  s-1. 
In  contrast,  ebb  tide  profiles  had  lower  diffusivities  at  mid-depth  and  higher  diffusivities 
near  the  surface  and  bottom.  These  asymmetries  were  captured  well  by  the  smooth  surfaces 
Ks  that  were  fitted  to  Kp  and  used  in  the  advection-diffusion  model. 
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4.4.2  Plankton  Profiles 


The  10-hour  series  of  plankton  samples  showed  that  larvae  were  temporally  patchy  even 
over  short  time  scales.  When  the  size  classes  were  pooled,  the  depth-averaged  larval  con¬ 
centration  (#Z/-1)  usually  varied  by  an  order  of  magnitude  or  more  within  each  10-hour 
period.  This  was  true  for  major  taxa  (Fig.  4-8)  and  at  the  genus  level  for  gastropods  (Fig. 
4-9).  Patches  of  Crepidula  spp.  and  Anachis  spp.  appeared  to  coincide,  whereas  Hycinassa 
were  most  abundant  when  the  other  2  genera  were  scarce  or  absent  (Fig.  4-9). 

Larvae  were  rarely  distributed  evenly  in  vertical  profiles  (Figs.  4-10  to  4-14).  Two 
plankton  profiles  were  omitted  because  they  were  collected  during  periods  of  slight  strat¬ 
ification  (22  July,  16:30  and  29  July,  10:30).  For  the  remaining  profiles  the  water  column 
was  very  well-mixed,  and  the  uneven  vertical  distributions  of  larvae  can  be  attributed  to 
behavior.  It  would  be  difficult,  however,  to  interpret  these  vertical  profiles  in  terms  of  be¬ 
havior  without  a  statistical  analysis.  Larval  vertical  distributions  depend  on  vertical  mixing 
profiles,  behavioral  responses  to  dissipation  rate,  and  the  dissipation  rate  profiles  as  well 
as  the  previous  vertical  distributions  (before  samples  were  collected).  The  maximum  like¬ 
lihood  analysis  (results  below)  allowed  us  to  extract  quantitative  behavior  estimates  from 
observed  larval  distributions  and  turbulence  profiles. 

4.4.3  Maximum-Likelihood  Behavior  Estimates 

The  response  of  crab  larvae  to  turbulence  was  significantly  different  (a  =  0.05)  on  flood  and 
ebb  tides  (Fig.  4-15a-b).  On  flood  tides,  larvae  had  a  nearly  constant  downward  velocity 
of  about  w  =  —0.1  cm  s-1.  On  ebb  tides,  larval  velocity  was  positively  related  to  the 
dissipation  rate,  so  that  larvae  were  swimming  upwards  at  e  greater  than  '--0.1  cm2  s“3. 

The  behavior  of  pooled  gastropod  larvae  was  also  significantly  different  on  flood  and  ebb 
tides.  For  the  large  size  fraction,  sample  sizes  were  very  small  and  b\  had  wide  confidence 
intervals.  As  a  result,  behavioral  estimates  for  the  two  size  fractions  were  never  significantly 
different.  Results  are  shown  for  both  size  fractions  combined  in  Fig.  4-15,  with  all  MLE’s 
given  in  Table  4.1.  For  pooled  gastropods,  larval  velocity  w  had  a  near-zero  slope  ( b\  = 
0.0125)  vs.  e  on  flood  tides  and  a  more  positive  slope  (b\  =  0.075)  vs.  £  on  ebb  tides  (Fig. 
4-15c-d).  The  estimates  of  w  were  below  the  neutral  buoyancy  line  for  almost  all  conditions. 

Gastropod  larvae  showed  genus-level  differences  in  the  response  to  turbulence.  Ilyanassa 
larvae  showed  almost  no  relationship  between  w  and  £,  with  negative  velocities  under  all 
conditions  (flood  tides  bo  =  —0.14  cm  s_1,  ebb  tides  6q  =  —0.16  cm  s_1).  There  was  a 
strong  negative  relationship  between  w  and  £  for  the  large  size  fraction  on  ebb  tides,  but 


66 


Table  4.1:  Maximum-likelihood  parameter  estimates  for  the  response  of  larvae  to  turbulence, 
w  =  bo  H~  b\  log^Q  £. 


taxon 

tidal  stage 

small 

fraction 

b0  b\ 

—  1  2  —1 
cm  s  s  cm  1 

large 

fraction 

bo  fci 

—  1  2  —1 

cm  s  1  s  cm  1 

both 

fractions 

60  Si 

—  1  2  —1 

cm  s  1  cm  A 

crabs 

flood 

-0.14 

-0.075 

-0.1 

0.0125 

-0.12 

0.0 

ebb 

0.2 

0.2 

0.08 

0.125 

0.12 

0.15 

gastropods 

flood 

-0.08 

0.0125 

-0.14 

0.025 

-0.08 

0.0125 

ebb 

-0.02 

0.075 

-0.02 

0.1 

-0.02 

0.075 

Ilyanassa 

flood 

-0.12 

0.0125 

-0.14 

-0.0125 

-0.14 

0.0125 

ebb 

-0.16 

-0.0125 

-0.3 

-0.225 

-0.16 

-0.0125 

Crepidula 

flood 

0.2 

0.125 

0.2 

0.125 

ebb 

0.06 

0.1 

0.2 

0.25 

0.06 

0.1 

Anachis 

flood 

-0.3 

-0.075 

0.2 

0.175 

-0.22 

-0.0375 

ebb 

-0.02 

0.1 

0.2 

0.2 

0.02 

0.125 
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large  confidence  interval  on  the  slope  b\  included  both  positive  and  negative  values  (Fig. 
4-16). 

riyanassa  larvae  in  the  large  size  fraction  were  likely  to  be  competent  for  metamor 
phosis,  and  were  expected  to  behave  similarly  to  competent  mud  snail  larvae  in  laboratory 
turbulence.  Field  estimates  of  II yanassa  behavior  in  turbulence  are  somewhat  different  than 
laboratory  results  (Fig.  4-16).  Most  of  the  laboratory  data  points  are  outside  the  confidence 
intervals  for  b\  in  both  flood  and  ebb  tide  estimates.  Maximum  likelihood  estimates  of  w 
are  below  the  laboratory  values  for  all  but  the  highest  turbulence  levels  ( e  >  10“03  cm2 
s-3). 

Both  Crepidula  spp.  and  Anachis  spp.  showed  a  positive  response  to  turbulence  on 
ebb  tides,  swimming  upward  in  strong  turbulence.  These  larvae  had  strong  responses  to 
turbulence  on  flood  tide  too,  but  the  confidence  intervals  for  61  included  both  positive  and 
negative  values. 

The  goodness  of  fit  tests  indicated  that  the  linear  behavior  model  (4.6)  was  incorrect  for 
all  larval  groups.  The  linear  model  could  not  be  rejected  for  large  Crepidula  larvae,  but  the 
sample  sizes  were  too  small  for  the  test  to  be  valid.  Although  larvae  had  a  strong  response 
to  turbulence,  the  shape  of  the  true  response  curve  is  probably  nonlinear  and  may  depend 
on  detection  limits  as  well  as  viscosity  and  other  physical  factors. 


4.5  Discussion 

4.5.1  Hydrodynamics 

Although  Barnstable  Harbor  is  a  sheltered  inlet,  it  is  very  turbulent,  probably  due  to  bottom 
friction  over  large-  and  small-scale  roughness  elements.  The  near-bottom  dissipation  rates 
of  up  to  e  =  101  cm2  s-3  are  higher  than  normal  for  tidal  channels  [14,  23,  43],  and  within 
the  range  of  dissipation  rates  in  the  surf  zone  [13,  42].  Near-bottom  drag  coefficients  were 
also  high  at  C&  —  0.01.  C&  is  typically  about  3  x  10~3  for  flow  over  smooth  sandy  substrates 
[15,  18],  but  the  present  estimates  of  C d  are  on  the  order  of  drag  coefficients  reported  over 
mobile  bed  forms  in  a  tidal  inlet  [46]  and  in  the  surf  zone  [31].  The  high  estimates  of  C d  in 
Barnstable  Harbor  are  probably  due  to  form  drag  over  large  dunes  (~  2  m  height,  ~  6  to 
10  m  wavelength).  Additional  bottom  roughness  features  include  smaller  sand  ripples  (~  1 
cm)  and  dense  aggregations  of  mud  snails  (~  1  cm).  Near-bottom  turbulence  around  these 
multi-scale  roughness  elements  is  probably  an  important  issue  for  larval  settlement,  but  is 
beyond  the  scope  of  this  paper. 
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Flows  were  tidally  asymmetric,  but  whereas  near-bottom  current  velocities  were  higher 
on  flood  tides,  u*,  C and  ew  were  higher  on  ebb  tides.  Tidal  asymmetries  in  turbulence 
were  damped  with  increasing  height  above  the  bottom,  and  could  reflect  the  fact  that  form 
drag  is  greater  on  ebb  tides  when  the  current  direction  opposes  the  direction  of  the  sand 
waves.  Larvae  are  expected  to  attach  to  the  bottom  more  successfully  when  shear  velocity 
is  low,  and  the  fact  that  near-bottom  shear  velocities  were  higher  on  ebb  tides  suggests  that 
flood  tides  would  be  a  more  opportune  time  for  larval  settlement. 

Turbulence  in  the  harbor  was  assumed  to  be  laterally  homogenous  for  the  analysis  of 
plankton  distributions,  but  this  assumption  could  be  wrong  given  the  presence  of  large  bed 
forms.  Sandwaves  probably  contribute  to  spatial  variability  in  turbulence  profiles  over  at 
least  the  wavelength  of  the  sandwaves  (~  6  to  10  m).  The  current  meters  were  located 
on  the  ridge  of  a  sandwave,  but  the  bottom  conditions  were  unknown  at  the  plankton 
mooring.  Therefore  there  is  some  uncertainty  in  vertical  turbulence  profiles  used  in  the 
plankton  analysis.  It  is  difficult  to  asses  the  potential  errors  in  behavior  estimates  due  to 
lateral  inhomogeneity  of  turbulence,  without  a  detailed  sensitivity  analysis.  Both  sites  were 
located  in  a  straight  channel  segment  though,  and  the  diffusivity  and  dissipation  estimates 
used  in  the  plankton  analysis  are  probably  correct  to  order  of  magnitude. 

4.5.2  Plankton  Variability 

Larval  concentrations  were  temporally  variable  over  10  h  sampling  periods,  and  this  vari¬ 
ability  is  attributed  to  the  horizontal  advection  of  larval  patches  past  the  sampling  site. 
The  maximum  spatial  extent  of  these  patches  can  be  estimated  by  integrating  the  surface 
current  velocity  over  each  sampling  period.  In  this  way,  22  July  samples  represent  4.0  km, 
27  July  samples  represent  4.6  km,  and  29  July  samples  represent  6.5  km.  These  estimates 
are  on  the  order  of  patch  sizes  reported  for  blue  crab  larvae  offshore  [25]  and  fiddler  crab 
larvae  downstream  of  a  spawning  site  [28].  Patches  that  were  present  for  shorter  time  peri¬ 
ods  (e.g.  29  July  crabs,  Fig.  4-8c,  and  29  July  Ilyanassa ,  Fig.  4-9b)  were  almost  certainly 
less  extensive  than  the  maximum  spatial  estimates.  At  low  tide,  a  larva  would  have  to 
travel  about  4.6  km  from  the  plankton  mooring  to  reach  the  outlet  to  Cape  Cod  Bay  or 
vice  versa.  Samples  representing  this  spatial  scale  had  order-of-magnitude  variability  in  lar¬ 
val  concentration,  indicating  that  larvae  are  distributed  heterogeneously  across  the  harbor. 
The  horizontal  patchiness  of  planktonic  larvae  could  contribute  to  patchy  settlement. 

The  characteristics  of  the  plankton-sample  time-series  allow  some  conjectures  about 
whether  larvae  were  spawned  locally  or  dispersed  to  the  harbor  from  elsewhere.  Crab  larva 
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were  most  abundant  on  mid-ebb  tides,  and  were  mostly  small,  early  stage  zoea  that  probably 
originated  within  Barnstable  Harbor.  Ilyanassa  larvae  probably  also  originated  within  the 
harbor,  given  that  they  were  most  concentrated  during  late  ebb  tides  and  that  the  adult 
snails  are  extremely  abundant  at  the  study  site.  In  contrast  Anachis  spp.  and  Crepidula 
spp.  were  most  concentrated  on  late  flood  and  early  ebb  tides.  Adults  of  these  species 
are  more  commonly  found  on  sandy  and  gravelly  beaches,  and  these  larvae  probably  were 
spawned  outside  the  harbor. 

4.5.3  Larval  Behavior 

This  analysis  produced  estimates  of  larval  behavior  that  are  reasonable  for  all  physically 
realistic  values  of  the  dissipation  rate,  despite  the  fact  that  larval  velocity  is  unbounded 
by  the  linear  model  (4.6).  Larvae  from  coastal  populations  probably  encounter  dissipation 
rates  ranging  from  10-7  cm2  s-3  in  the  open  ocean  [6]  to  103  cm2  s~3  in  the  surf  zone 
[13].  Using  some  of  the  most  extreme  MLE’s  (bo  =  0.2  cm  s~!and  —  0.2  s2  cm’1  for 
small  crab  larvae  on  ebb  tides,  Table  4.1),  the  average  larval  velocity  would  span  -1.2  to 
0.8  cm  s-1  over  the  entire  range  of  potential  turbulence  conditions.  These  estimates  are 
well  within  the  range  of  swimming  and  sinking  velocities  for  crab  larvae  [reviewed  in  2] ,  and 
within  reported  limits  for  veliger  larvae  [e.g.,  9,  19]. 

Although  the  velocity  estimates  were  reasonable,  the  linear  response  to  turbulence  (4.6) 
was  an  incorrect  model  for  all  larval  groups.  This  result  is  not  surprising.  The  linear 
model  assumes  that  larvae  are  equally  able  to  detect  strong  and  weak  turbulence,  and  that 
larval  behavior  depends  only  on  turbulence.  Larval  detection  of  turbulence  almost  certainly 
depends  on  the  length  scale  rj  of  the  smallest  eddies  relative  to  larval  size.  The  smallest 
eddies  are  smaller  than  larvae  in  strong  turbulence  (for  e  —  101  cm2  s“3,  rj  ^  0.2  mm) 
and  larger  than  larvae  in  weak  turbulence  (for  e  —  10-4  cm2  s~3,  rj  ^  3  mm).  There  is 
probably  some  turbulence  threshold  that  is  below  larval  detection  limits.  It  is  also  very  likely 
that  larvae  respond  to  additional  environmental  cues  that  vary  non-linearly  with  turbulence 
(discussed  below).  Even  though  the  linear  model  is  incorrect,  the  behavior  estimates  confirm 
that  turbulence  is  an  important  cue  for  behavior  of  crabs  and  gastropods  in  the  field. 

Crab  larvae 

Behavior  estimates  for  crab  larvae  suggest  that  the  response  of  zoea  to  turbulence  enhances 
export  of  larvae  to  coastal  areas.  These  results  are  consistent  with  other  evidence  for 
selective  tidal  stream  transport  (STST)  of  crab  larvae.  The  distributions  of  Pachygrapsus 
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cmssipes  larvae  in  San  Diego  Bay  strongly  suggest  that  these  larvae  sink  on  flood  tides  and 
swim  upward  on  ebb  tides,  resulting  in  the  export  of  larvae  from  the  bay  [5].  Likewise  in  this 
study,  crab  larvae  swam  upwards  in  strong  turbulence  during  ebb  tides  and  had  a  constant 
downward  velocity  during  flood  tides.  The  net  result  would  be  the  export  of  larvae  from 
the  harbor  to  Cape  Cod  Bay.  Tidally-dependent  behaviors  are  species-specific,  and  not  all 
species  have  larval  behaviors  that  enhance  export  from  natal  estuaries  [3,  5].  In  this  study 
the  behavior  of  crab  larvae  was  analyzed  at  a  high  taxonomic  level,  yet  zoea  still  showed  a 
strong  export-enhancing  response  to  turbulence.  The  samples  could  have  been  dominated 
by  a  species  with  STST  behaviors.  Alternatively,  it  is  possible  that  the  most  common  crab 
species  in  the  harbor  share  similar  STST  behavioral  strategies. 

The  fact  that  crab  larvae  swim  up  on  ebb  tide  turbulence  but  not  on  flood  tide  turbulence 
indicates  that  behavior  is  determined  by  multiple  cues,  as  would  be  expected  for  STST. 
Blue  crab  ( Callinectes  sapidus)  megalopae  swim  upwards  in  turbulence,  but  the  response  is 
stronger  if  accompanied  by  an  increase  in  salinity  mimicking  flood  tide  [45].  In  the  present 
study  the  salinity  range  was  very  small  over  the  sampling  periods  (<—  1  psu  over  10  hours), 
but  temperature  changes  were  large  (~  6  °C  over  10  hours)  and  could  also  play  a  role  in 
tidally-dependent  behavior. 


Gastropod  larvae 

The  turbulence  response  of  gastropod  larvae  was  clearly  species-dependent,  and  could  in¬ 
dicate  that  larvae  from  different  habitat  types  use  turbulence  differently.  Ryanassa  larvae 
were  sinking  under  most  conditions.  This  behavior  would  enhance  larval  retention  in  the 
harbor  and  nearshore  zone  [30],  and  could  allow  larvae  to  stay  near  natal  habitats  through¬ 
out  the  development  period.  Large  Ryanassa  larvae  had  greater  sinking  velocities  on  strong 
ebb-tide  turbulence,  which  would  retard  the  export  of  potentially  competent  larvae  to  Cape 
Cod  Bay.  Large  Ryanassa  larvae  were  present  in  few  plankton  profiles  though,  and  the 
behavior  estimates  are  inconclusive  based  on  the  confidence  intervals  for  bj .  Based  on  these 
few  plankton  profiles,  mud  snail  larvae  responded  differently  in  the  field  than  they  did  in 
laboratory  turbulence,  indicating  that  behavior  depends  on  additional  cues. 

Crepidula  and  Anachis  larvae  behaved  more  like  crab  larvae,  swimming  upward  on 
strong  ebb  tide  turbulence.  Behaviors  were  inconsistent  on  flood  tides;  Crepidula  swam 
upward  and  Anachis  sank  on  flood  tide  turbulence.  Crepidula  larvae  apparently  avoid  the 
bottom  in  any  turbulent  conditions  whereas  Anachis  behavior  would  favor  export  from 
estuaries.  Adult  snails  of  these  species  inhabit  sandy  to  gravelly  beaches,  rather  than 
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intertidal  mud  flats.  The  turbulence  responses  of  these  larvae  could  allow  them  to  avoid 
setthng  in  unfavorable  estuarine  habitats. 

One  unexpected  result  of  the  behavior  estimates  was  that  larvae  were  sinking  under  a 
broad  range  of  conditions,  even  when  there  was  no  response  to  turbulence  (Fig.  4-15a,c,e). 
This  result  has  several  possible  explanations.  As  mentioned  above,  sinking  could  allow 
larval  retention  near  favorable  habitats.  Sinking  could  also  be  a  response  to  predators,  as 
discussed  in  Ch.  2,  but  it’s  equally  likely  that  sinking  was  due  to  negative  phototaxis.  All 
of  the  samples  were  collected  on  bright  sunny  days  in  shallow  water,  and  widespread  sinking 
could  be  a  result  of  surface-avoidance. 

4.5.4  Summary 

On  average,  turbulence  was  an  important  behavioral  cue  for  all  larvae  sampled,  at  least 
dining  ebb  tides.  Larvae  of  crabs  and  beach-dwelling  gastropods  swam  upward  on  ebb¬ 
tide  turbulence,  a  behavior  that  would  enhance  larval  transport  out  of  the  harbor.  In 
contrast  large  mud  snail  larvae  sank  on  ebb-tide  turbulence,  a  behavior  that  would  promote 
larval  retention.  No  turbulence  response  was  detectable  for  pooled  gastropods  and  crab 
larvae  on  flood  tides.  This  negative  result  could  indicate  the  presence  of  multiple  species 
whose  turbulence  responses  cancelled  out,  or  the  suppression  of  a  turbulence  response  by 
additional  cues  (e.g.  increasing  salinity  [45]).  These  results  support  the  hypothesis  that 
larval  behavior  in  turbulence  contributes  to  the  supply  of  larvae  to  benthic  habitats.  Genus- 
specific.  responses  to  turbulence  suggest  that  larvae  could  select  habitats  over  relatively  large 
spatial  scales  (>  km). 
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Figure  4-3:  (a-c)  Temperature  and  (d-f)  salinity  plots  for  plankton  sampling  periods  on 
(a,  d)  22  July,  (b,  e)  27  July,  and  (c,  f)  29  July.  Double  lines  represent  down*  and  up-casts. 
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Figure  4-4:  July  measurements  of  (a)  surface  current  velocity  (positive  on  ebb  tide,  neg¬ 
ative  on  flood  tide),  (b)  depth-averaged  shear  velocity,  (c)  depth- averaged  drag  coefficient, 
and  (d)  surface  dissipation  rate  from  ADCP  measurements  (z  =  1.5  m  to  surface). 
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Figure  4-5:  July  measurements  of  (a)  current  velocity  (positive  on  ebb  tide,  negative 
on  flood  tide),  (b)  shear  velocity,  (c)  drag  coefficient,  and  (d)  dissipation  rate  from  ADV 
measurements  {z  =  0.7  m).  Missing  data  points  were  discarded  because  of  excessive  noise 
due  to  sensor  obstruction. 
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Figure  4-6:  (a-c)  Measured  eddy  diffusivity  Kp  and  (d-f)  fitted  eddy  diffusivity  K„  for 
plankton-sampling  periods  of  (a,  d)  22  July,  (b,  e)  27  July,  and  (c,  f)  29  July. 
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Figure  4-7:  (a-c)  Measured  dissipation  rate  ep  and  (d-f)  fitted  dissipation  rate  e3  for 
plankton-sampling  periods  of  (a,  d)  22  July,  (b,  e)  27  July,  and  (c,  f)  29  July. 
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Figure  4-8:  Time  series  of  average  total  larval  concentration  (#L-1)  per  profile  for  (a)  22 
July,  (b)  27  July,  and  (c)  29  July.  Size  classes  are  pooled  and  larvae  are  grouped  by  major 
taxa.  Asterisks  indicate  excluded  profiles.  Background  shading  indicates  ebb  tides.  Error 
bars  are  1  standard  error. 
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Figure  4-9:  Time  series  of  average  total  larval  concentration  (#L-1)  per  profile  of  gas¬ 
tropods  for  (a)  27  July,  and  (b)  29  July  (22  July  gastropods  were  not  identified  to  species). 
Size  classes  are  pooled  and  larvae  are  grouped  by  genera.  Asterisks  indicate  excluded  pro¬ 
files.  Background  shading  indicates  ebb  tides.  Error  bars  are  1  standard  error. 
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Figure  4-10:  Ten-hour  series  of  vertical  plankton  profiles  of  crab  zoea  in  (o)  large  and  (•) 
small  size  fractions  for  (a)  22  July,  (b)  27  July,  and  (c)  29  July.  Larval  concentrations  are 
shown  as  proportion  of  the  profile  total,  which  is  indicated  above  each  profile.  Profiles  with 
fewer  larvae  than  samples  are  not  plotted. 
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Figure  4-11:  Ten-hour  series  of  vertical  plankton  profiles  of  pooled  gastropod  larvae  in 
(o)  large  and  (•)  small  size  fractions  for  (a)  22  July,  (b)  27  July,  and  (c)  29  July.  Larval 
concentrations  are  shown  as  proportion  of  the  profile  total,  which  is  indicated  above  each 
profile.  Profiles  with  fewer  larvae  than  samples  are  not  plotted. 
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Figure  4-12:  Ten-hour  series  of  vertical  plankton  profiles  of  Ilyanassa  obsoleta  larvae  in 
(o)  large  and  (•)  small  size  fractions  for  (a)  27  July,  and  (b)  29  July.  Larval  concentrations 
are  shown  as  proportion  of  the  profile  total,  which  is  indicated  above  each  profile.  Profiles 
with  fewer  larvae  than  samples  are  not  plotted. 
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Figure  4-13:  Ten-hour  series  of  vertical  plankton  profiles  of  Crepidula  spp.  larvae  in  (o) 
large  and  (•)  small  size  fractions  for  (a)  27  July,  and  (b)  29  July.  Larval  concentrations  are 
shown  as  proportion  of  the  profile  total,  which  is  indicated  above  each  profile.  Profiles  with 
fewer  larvae  than  samples  are  not  plotted. 
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Figure  4-14:  Ten- hour  series  of  vertical  plankton  profiles  of  Anachis  spp.  larvae  in  (o) 
large  and  (•)  small  size  fractions  for  (a)  27  July,  and  (b)  29  July.  Larval  concentrations  are 
shown  as  proportion  of  the  profile  total,  which  is  indicated  above  each  profile.  Profiles  with 
fewer  larvae  than  samples  are  not  plotted. 
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Figure  4-15:  Maximum  likelihood  estimates  of  larval  velocity  w  =  bo  +  &ilog10£  vs. 
dissipation  rate  e  (solid  lines).  Small  and  large  size  classes  are  combined.  Separate  estimates 
are  shown  for  (a,c,e,g,i)  flood  tide  and  (b,d,f,h,j)  ebb  tide  samples.  Dotted  lines  indicate 
neutral  buoyancy,  dashed  lines  indicate  a  95%  confidence  interval  for  the  slope  &i. 
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Figure  4-16:  Maximum  likelihood  estimates  of  larval  velocity  w  =  feo  +  bilog10£  vs- 
dissipation  rate  £  for  large  llyanassa  larvae,  with  laboratory  data  superimposed  (redrawn 
from  [9]).  Separate  estimates  are  shown  for  (a)  flood  tide  and  (b)  ebb  tide  samples.  Dotted 
lines  indicate  neutral  buoyancy,  dashed  lines  indicate  a  95%  confidence  interval  for  the  slope 

h. 
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Chapter  5 


Summary  and  Conclusions 


5.1  Summary  of  Results 

The  goals  of  this  thesis  were  to  quantify  the  behavior  of  gastropod  larvae  (mud  snails 
llyanassa  obsoleta)  in  turbulence,  and  to  investigate  how  that  behavior  affects  larval  supply 
in  a  turbulent  coastal  inlet.  I  used  several  approaches,  including  laboratory,  modeling,  and 
field  studies,  to  explore  this  problem.  The  laboratory  study  was  necessary  to  quantitatively 
describe  the  response  of  larvae  to  turbulence.  The  model  provided  insights  on  the  theoretical 
effects  of  this  behavior  on  larval  supply  and  settlement.  The  field  study  produced  some 
discoveries  about  larval  responses  to  turbulence  that  would  have  been  very  difficult  to 
obtain  in  the  laboratory.  My  results  suggest  that  larval  responses  to  turbulence  significantly 
affect  larval  supply  to  benthic  habitats,  and  that  larvae  from  multiple  habitats  could  use 
turbulence  as  a  cue  for  habitat  selection. 

In  Chapter  2,  I  quantified  the  behavior  of  competent  mud  snail  larvae  in  laboratory 
turbulence.  The  larvae  had  three  behavioral  modes:  swimming,  hovering,  and  sinking. 
The  proportion  of  larvae  in  the  sinking  mode  increased  exponentially  with  the  turbulence 
dissipation  rate.  As  a  result,  the  average  larval  velocity  shifted  from  positive  (upward)  to 
negative  (downward)  at  a  dissipation  rate  of  about  e  ~  10-1  cm2  s~3.  This  turbulence 
threshold  occurs  when  the  Kolmogorov-scale  eddies  are  about  the  size  of  the  larvae,  sup¬ 
porting  the  notion  that  larvae  detect  turbulence  by  interacting  with  the  smallest  eddies. 
The  shift  in  behavior  occurred  at  a  turbulence  level  found  only  in  coastal  areas,  and  the 
response  to  turbulence  is  expected  to  influence  the  supply  of  larvae  to  benthic  habitats. 

In  Chapter  3,  I  used  a  model  to  characterize  how  turbulence-induced  sinking  affects 
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larval  supply  and  settlement  in  tidal  channels  or  inlets.  Behaving  larvae  (those  that  sink 
in  turbulence)  had  higher  larval  supply  than  passive  larvae  (those  with  constant  buoyancy) 
during  flood  and  ebb  tides,  but  not  during  slack  tides.  Higher  larval  supply  gave  behaving 
larvae  a  settlement  advantage  over  passive  larvae  under  some  hydrodynamic  conditions. 
In  strong  currents,  behaving  larvae  settled  more  successfully  than  passive  larvae  if  the 
settlement  velocity  was  non-zero  during  flood/ebb  tides.  Behaving  larvae  settled  more 
successfully  in  stronger  currents  than  in  weaker  ones.  These  results  suggest  that  sinking  in 
turbulence  makes  larvae  more  likely  to  settle  in  turbulent  inlets  such  as  Barnstable  Harbor. 

In  Chapter  4,  I  used  data  from  Barnstable  Harbor  to  estimate  how  larvae  respond 
to  turbulence  in  the  field.  When  larval  velocity  was  estimated  as  a  linear  function  of  the 
dissipation  rate,  larvae  showed  a  strong  response  to  turbulence.  For  all  larval  groups  (crabs, 
gastropods,  mud  snails,  slipper  shells,  and  dove  shells)  behavior  differed  on  flood  and  ebb 
tides,  indicating  that  other  environmental  cues  are  also  important  determinants  of  behavior. 
Gastropod  larvae  had  genus-specific  behaviors,  suggesting  that  turbulence  provides  a  cue 
for  large-scale  habitat  selection  by  species  from  multiple  habitat  types. 

5.2  Evolutionary  Context 

There  are  two  main  hypotheses  for  why  larvae  sink  in  turbulence.  Until  recently  it  was 
generally  assumed  that  larvae  sank  to  avoid  predators  [e.g.,  9].  In  this  thesis  I  propose  that 
larvae  of  coastal  species  use  turbulence  as  an  indicator  of  potentially  favorable  habitats,  and 
that  sinking  allows  larvae  to  get  to  the  bottom  and  explore  for  settlement  sites  in  coastal 
areas.  The  predator-avoidance  hypothesis  and  settlement  hypothesis  are  not  contradictory, 
making  it  difficult  to  test  or  reject  either  hypothesis. 

Sinking  in  turbulence  probably  increases  larval  fitness  in  multiple  ways.  Larvae  that 
sink  in  turbulence  may  escape  being  eaten  by  predators  [but  see  4] ,  and  could  have  lower 
mortality  rates  during  dispersal.  As  shown  in  Chapter  3,  larvae  that  sink  in  turbulence 
are  also  more  likely  to  settle  into  favorable  habitats  than  into  unfavorable  ones,  and  should 
have  lower  mortality  rates  at  or  after  settlement.  The  combined  reductions  in  mortality 
rate  probably  reinforce  the  adaptive  benefits  of  sinking  in  turbulence. 

Some  larvae  may  sink  “mistakenly,”  for  example  in  storm-generated  turbulence  at  sea, 
and  it  is  unknown  whether  the  costs  outweigh  the  benefits  for  larvae  that  sink  in  storms. 
Mud  snail  larvae  sink  at  a  turbulence  threshold  that  is  higher  than  the  dissipation  rates 
measured  during  mild  offshore  storms  (up  to  e  =  10~2  cm2  s~3  in  winds  up  to  16  ra  s'"1 
[3,  5]).  Dissipation  rates  in  stronger  storms  would  probably  reach  levels  well  above  the 
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threshold  for  larval  sinking.  It  is  expected  that  larvae  would  sink  to  calmer  depths  during 
storms  [e.g.,  1],  but  this  sinking  might  be  costly  if  sinking  slows  the  return  transport  of  lar¬ 
vae  to  coastal  habitats  and  prevents  larvae  from  settling.  Although  sinking  in  storms  might 
be  counterproductive  for  settlement,  it  could  be  beneficial  for  avoiding  predators.  Storm¬ 
generated  turbulence  would  increase  contact  rates  between  larvae  and  predators  [e.g.,  6,  7], 
and  sinking  would  still  allow  larvae  to  reach  calmer  depths  where  contact  rates  with  preda¬ 
tors  are  lower.  Any  costs  of  sinking  in  storm-related  turbulence  are  probably  outweighed 
overall  by  the  benefits  of  avoiding  predators  and  settling  more  successfully. 


5.3  Unanswered  Questions  and  Future  Work 

5.3.1  Laboratory  vs.  field  estimates  of  behavior 

Responses  of  mud  snail  larvae  to  turbulence  were  estimated  both  in  the  laboratory  and 
in  the  field.  In  the  laboratory,  average  larval  vertical  velocities  switched  from  upward  to 
downward  at  e  =  10-1  cm2  s~3.  In  the  field,  larvae  responded  differently  to  turbulence  on 
ebb  and  flood  tides,  and  only  the  ebb-tide  behavior  estimates  were  similar  to  laboratory 
estimates.  For  both  tidal  stages,  field  estimates  of  larval  velocity  were  lower  than  laboratory 
estimates  except  at  the  highest  turbulence  levels  [e  >  1  cm2  s~3). 

There  are  several  potential  reasons  for  differences  between  field  and  laboratory  estimates 
of  larval  velocity:  1.  Other  physical  factors  (e.g.  salinity,  temperature)  besides  turbulence 
varied  in  the  field,  and  larval  behavior  could  depend  on  additional  cues.  2.  Temperature 
changes  in  the  field  could  have  affected  larval  swimming  abilities  by  altering  viscosity  and 
by  altering  larval  biochemical  rates.  3.  Laboratory  estimates  of  larval  velocity  decreased 
exponentially  with  turbulence;  the  linear  response  model  used  in  the  field  study  was  inca¬ 
pable  of  capturing  non-linear  responses  to  turbulence.  4.  Field  data  for  large  mud  snail 
larvae  were  sparse,  resulting  in  wide  confidence  intervals  on  behavior  estimates.  The  dif¬ 
ferences  between  laboratory  and  field  estimates  could  be  due  partly  to  differences  in  larval 
behavior  under  artificial  vs.  natural  conditions,  and  due  partly  to  the  difficulty  in  obtaining 
good  behavior  estimates  from  field  data.  This  uncertainty  could  be  resolved  partially  by 
using  more  plankton  data  and  by  fitting  non-linear  behavior  models  to  the  field  plankton 
distributions. 
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5.3.2  Omissions  of  the  model 


The  model  in  Chapter  3  is  a  simplification  of  bound  ary-  layer  hydrodynamics  and  of  larval 
behavior.  This  simplified  model  provides  insights  into  the  effects  of  turbulence-induced 
sinking  on  larval  supply  and  settlement,  but  many  questions  remain  about  larval  settlement 
in  the  real  world.  For  example,  mud  snail  larvae  have  sediment-selection  behaviors  in 
still  water  [8]  -  are  these  behaviors  effective  in  turbulent  tidal  inlets?  If  so  what  are  the 
combined  effects  of  turbulence-induced  sinking  and  sediment  selection  behavior  on  larval 
settlement?  In  the  field,  larval  behavior  differs  on  flood  and  ebb  tides  -  how  does  tidal- 
stage  dependent  behavior  affect  temporal  patterns  of  larval  supply?  Other  questions  relate 
to  the  characteristics  of  the  tidal  boundary  layer.  What  happens  to  larvae  that  approach 
the  bottom  and  become  entrained  in  mobile  sediments?  What  about  larvae  settling  at  the 
water’s  edge?  How  do  turbulence  intermittency  and  calm  eddies  around  roughness  elements 
affect  the  ability  of  larvae  to  settle  during  peak  flood  and  ebb  tides?  These  questions  were 
beyond  the  scope  of  this  thesis  but  could  be  answered  with  additional  laboratory,  field,  and 
modeling  studies. 


5.3.3  Future  work 

Many  questions  remain  about  the  ecological  benefits  of  larval  sinking  in  turbulence,  and 
some  of  these  questions  could  be  answered  with  future  modeling  efforts.  For  example,  how 
does  the  response  to  turbulence  affect  settlement  success  in  patchy  environments?  And  do 
larvae  benefit  more  from  settlement  responses  to  cues  in  the  water  column  or  very  near 
the  bottom  (i.e.  substrate-selection  behavior)?  I  hope  to  address  these  questions  with  a 
2-dimensional  model  that  incorporates  substrate  patchiness  and  larval  substrate-selection 
behavior. 

Thus  far,  horizontal  larval  dispersal  has  been  treated  generally  as  a  passive  process. 
Larval  behavior  in  turbulence  appears  to  influence  vertical  distributions  and  large-scale  (> 
km)  settlement  patterns,  and  probably  also  affects  patterns  of  dispersal.  It  is  now  possible 
to  predict  patterns  of  passive  larval  dispersal  using  sophisticated  physical  oceanographic 
models  (e.g.  FVCOM  [2]).  However,  physical  transport  models  typically  ignore  larval 
behavior  and  contain  large  biological  uncertainties.  Uncertainties  in  dispersal  estimates 
could  be  reduced  by  incorporating  quantitative  descriptions  of  larval  behavior  in  physical 
transport  models. 
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5.4  Significance 

In  this  thesis  I  use  two  new  approaches  for  estimating  larval  behaviors  in  turbulence.  The 
mixture  model  (Chapter  2)  makes  it  possible  to  untangle  larval  behavior  from  fluid  motion 
in  laboratory  flows.  With  some  further  development  of  the  model  and  parameter-estimation 
algorithm,  this  approach  could  be  used  to  estimate  plankton  behaviors  from  field-collected 
data  including  video  archives  from  marine  observatories.  The  maximum-likelihood  analysis 
of  larval  distributions  in  the  field  (Chapter  4)  provides  another  means  of  untangling  larval 
behavior  from  turbulent  flows.  Drawbacks  to  this  approach  are  the  sparseness  of  biological 
data  and  the  necessary  separation  distance  between  sites  of  biological  and  physical  data 
collection.  Despite  these  difficulties,  this  method  of  analysis  provides  quantitative  behavior 
estimates  that  are  an  improvement  over  the  more  speculative  inferences  of  the  past. 

This  research  supports  the  hypothesis  that  larvae  use  turbulence  as  a  cue  to  sink  to  the 
bottom  in  potentially-favorable  coastal  habitats.  It  is  still  impossible  to  reject  the  hypothesis 
that  larvae  sink  to  avoid  predators,  but  there  is  no  reason  to  believe  that  the  two  hypotheses 
are  contradictory.  Modeling  results  (Chapter  3)  show  clearly  that  sinking  in  turbulence  will 
increase  the  supply  of  larvae  to  benthic  habitats  in  turbulent  tidal  inlets.  Larvae  that  sink  in 
turbulence  are  more  likely  to  settle  in  favorable  habitats  such  as  Barnstable  Harbor  than  in 
unfavorable,  calmer  areas.  Genus-specific  responses  to  turbulence  (Chapter  4)  suggest  that 
the  use  of  turbulence  for  habitat  selection  could  be  widespread  among  coastal  gastropods. 
These  results  demonstrate  that  larval  behavior  is  important  and  should  no  longer  be  ignored 
in  models  of  dispersal  and  settlement. 
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S inking  behavior  of  gastropod  larvae  ( llyanassa  obsoleta)  in  turbulence 


Heidi  L.  Fuchs,  Lauren  S.  Mullineaux,  and  Andrew  R.  Solow 


Web  appendix  1.  Modifled  expectation-maximization 
algorithm 

For  behavioral  experiments  in  turbulence,  the  mixing  pro¬ 
portions  aT,  and  the  fluid  velocity  variance  of  were  esti¬ 
mated  by  the  expectation-maximization  (EM)  algorithm  be¬ 
low.  We  used  measured  values  for  the  mean  flow  velocity 
fia  and  maximum  likelihood  estimates  (from  still-water  data) 
for  the  behavioral  velocity  means  fi,  and  variances  of  of 
larvae  in  each  mode.  The  relative  mean  velocity  Mr'  for  each 
mode  was  known  (/ir.  —  fxa  +  fir).  Here,  i  =  (1  :  g)  is  the 
mixture  component  index  for  g  =  3  modes,  j  =  (1  :n)  is  the 
data  value  index  for  n  measurements,  and  k  is  the  iteration 
index.  The  probability  density  of  a  measured  velocity  wLj, 
given  that  it  was  drawn  from  a  norma)  distribution  with  pa¬ 
rameters  Or1  —  (Mr't  erf),  is  f(wVJ\  $$■*).  The  probability  that 
data  value  j  came  from  component  i,  given  the  parameters 
at  step  k ;  is  t,*.  The  expected  value  of  the  log-likelihood 
function,  given  the  parameters  at  step  k,  is  L{0T\  0f).  The 
algorithm  was  repeated  until  the  value  of  L  increased  by 
<10-4  per  iteration.  Only  Eqs.  4  and  5  differ  from  the  stan¬ 
dard  algorithm  described  in  McLachlan  and  Peel  (2000). 
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McLachlan,  G.,  and  D.  Peel.  2000.  Finite  mixture  models.  John 
Wiley  &  Sons. 
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Appendix  B 


Diffusivity  and  Dissipation 
Estimates  for  Maximum  Likelihood 
Analysis  (Chapter  4) 


For  maximum  likelihood  analysis  of  plankton  profiles,  the  diffusivity  and  dissipation  terms 
in  (4.5,  4.6)  had  to  be  accurate,  to  minimize  bias  in  the  estimates  of  larval  vertical  ve¬ 


locity,  as  well  as  continuous,  to  allow  numerical  solution  of  (4.5).  Simple  models  of  eddy 


diffusivity  and  dissipation  rate  were  poor  predictors  of  measured  values,  particularly  on  ebb 
tides.  Measured  diffusivity  and  dissipation  rates  were  calculated  at  10-minute  intervals  and 
were  discontinuous.  Therefore,  high-order  polynomial  functions  were  fitted  to  field  data  to 
produce  diffusivity  and  dissipation  functions  that  were  both  accurate  and  continuous  for 
use  in  (4.5). 

Smooth  diffusivity  functions  Ks  =  f(z/H,  t)  were  fitted  to  measured  diffusivity  Kp  for 
each  plankton  sampling  period.  No  near-bottom  diffusivity  measurements  were  available 
because  ADCP  measurement  bins  began  at  zt>  =  1.5  m.  In  order  to  fit  Ks  over  the  entire 
water  column,  the  diffusivity  was  extrapolated  from  —  1.5  down  to  z  —  0.  The  diffusivity 
between  the  lowest  bin  and  the  bottom  was  estimated  in  25  cm  increments  by 


(B.l) 


where  the  Reynolds  stress  was  extrapolated  to  the  bottom  using 


(B.2) 
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The  along-channel  velocity  shear  was  extrapolated  to  the  bottom  by  assuming  that  U (z) 
decreases  linearly  vs.  In  £  to  t/(0)  =  0  (the  no-slip  condition),  so 


dU*  _  U(zb)  1 
dz  In  z\y  z 


(B.3) 


A  polynomial  response  surface  was  fitted  to  Kv(z/H,t)  on  a  logio  scale  to  produce  a 
continuous-time  diffusivity  function  Ks.  The  polynomial  orders  of  2  and  t  were  selected 
as  those  that  provided  the  best  fit  to  the  data  for  each  sampling  period  (Table  B.l). 

Smoothed  field  data  were  used  likewise  for  the  turbulence  dissipation  rate  in  (4.6).  The 
dissipation  rate  between  the  lowest  bin  (zb  =  1.5  m)  and  the  bottom  was  estimated  in  25 
cm  increments  by 


£p  ~ 


KZ 


-±) 

HJ 


(B.4) 


where  k  is  von  Karmann’s  constant  (=  0.4).  A  response  surface  e3  =  f(z/Hrt)  was  fitted  to 
the  measured  dissipation  rate  ev  on  a  logio  scale  for  each  plankton  sampling  period,  using 
the  orders  of  z  and  t  that  provided  the  best  fit  (Table  B.l). 


Table  B.l:  Orders  of  z  and  t  that  provided  the  best  fits  of  Ks  and  es  to  measured  diffusivity 
Kp  and  dissipation  rate  ev  for  plankton  sampling  periods  (Ch.  4).  Also  given  are  R 2  for 
regressions  of  Ks  vs.  Kp  and  e3  vs.  ep. 


Ks 

sampling  period 

z  order 

t  order 

R2 

z  order 

t  order 

fF 

22  July 

20 

9 

0.62 

20 

7 

0.68 

27  July 

19 

6 

0.57 

20 

9 

0.70 

29  July 

20 

10 

0.59 

19 

6 

0.64 
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